Kawasaki dynamics with two types of particles: 

critical droplets 



F. den Hollander ^ ^ 
F.R. Nardi ^ 2 
A. Troiani ^ 

August 1, 2012 



Abstract 

This is the third in a series of three papers in which we study a two-dimensional lattice 
gas consisting of two types of particles subject to Kawasaki dynamics at low temperature 
in a large finite box with an open boundary. Each pair of particles occupying neighboring 
sites has a negative binding energy provided their types are different, while each particle has 
a positive activation energy that depends on its type. There is no binding energy between 
particles of the same type. At the boundary of the box particles are created and annihilated 
in a way that represents the presence of an infinite gas reservoir. We start the dynamics 
from the empty box and are interested in the transition time to the full box. This transition 
is triggered by a critical droplet appearing somewhere in the box. 

In the first paper we identified the parameter range for which the system is metastable, 
showed that the first entrance distribution on the set of critical droplets is uniform, computed 
the expected transition time up to and including a multiplicative factor of order one, and 
proved that the nucleation time divided by its expectation is exponentially distributed, all 
in the limit of low temperature. These results were proved under three hypotheses, and 
involved three model- dependent quantities: the energy, the shape and the number of critical 
droplets. In the second paper we proved the first and the second hypothesis and identified 
the energy of critical droplets. In the third paper we prove the third hypothesis and identify 
the shape and the number of critical droplets, thereby completing our analysis. 

Both the second and the third paper deal with understanding the geometric proper- 
ties of subcritical, critical and supercritical droplets, which are crucial in determining the 
metastable behavior of the system, as explained in the first paper. The geometry turns out 
to be considerably more complex than for Kawasaki dynamics with one type of particle, for 
which an extensive literature exists. The main motivation behind our work is to understand 
metastability of multi-type particle systems. 
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1 Introduction 



Motivation. The main motivation behind the present work is to understand mctastabiUty of 
multi-type particle systems subject to conservative stochastic dynamics. In the past ten years, a 
good understanding has been achieved of the metastable behavior of the lattice gas subject to 
Kawasaki dynamics, i.e., a conservative dynamics characterized by random hopping of particles 
of a single type with hardcore repulsion and nearest-neighbor attraction. The analysis was based 
on a combination of techniques from large deviation theory, potential theory, geometry and 
combinatorics. In particular, a precise description has been obtained of the time to nucleation 
(from the "gas phase" to the "liquid phase"), the shape of the critical droplet triggering the 
nucleation, and the typical nucleation path, i.e., the typical growing and shrinking of droplets. 
For an overview we refer the reader to two recent papers presented at the 12th Brazilian School of 
Probability: Gaudilliere and Scoppola [10] and Gaudilliere [9]. For an overview on metastability 
and droplet growth in a broader context, we refer the reader to the monograph by Olivieri and 
Vares [22], and the review papers by Bovier [3], [4], den Hollander [11], Olivieri and Scoppola [21]. 

The model we study constitutes a first attempt to generalize the results in Bovier, den Hollan- 
der and Nardi [6] for two-dimensional Kawasaki dynamics with one type of particle to multi-type 
particle systems. We take a large finite box A C Z^. Particles come in two types: type 1 and 
type 2. Particles hop around subject to hard-core repulsion, and are conserved inside A. At the 
boundary of A particles are created and annihilated as in a gas reservoir, where the two types 
of particles have different densities e~^^^ and e~^^^. We assume a binding energy U between 
particles of different type, and no binding energy between particles of the same type. Because 
of the "antiferromagnetic" nature of the interaction, configurations with minimal energy have a 
"checkerboard" structure. The phase diagram of this simple model is already very rich. The 
model can be seen as a conservative analogue of the Blume-Capel model investigated by Cirillo 
and Olivieri [8]. 

Our model describes the condensation of a low-temperature and low-density supersaturated 
lattice gas. We are interested in studying the nucleation towards the liquid phase represented 
by the checkerboard configuration ffl, starting from the gas phase represented by the empty 
configuration □. It turns out that the geometry of the energy landscape is much more complex 
than for the model of Kawasaki dynamics with one type of particle. Consequently, it is a 
somewhat delicate matter to capture the proper mechanisms behind the growing and shrinking 
of droplets. Our proofs use potential theory and rely on ideas developed in Bovier, den Hollander 
and Nardi [6] for Kawasaki dynamics with one type of particle. 

Two previous papers. In [14] we identified the values of the parameters for which the model 
properly describes the condensation of a supersaturated gas and exhibits a metastable behavior. 

Under three hypotheses, we determined the distribution and the expectation of the nucleation 
time, and identified the so-called critical configurations that satisfy a certain "gate property". The 
first hypothesis assumes that configuration ffl, corresponding to the liquid phase, is a minimizer 
of the Hamiltonian. The second hypothesis requires that the valleys of the energy landscape are 
not too deep. The third hypothesis requires that the critical configurations have an appropriate 
geometry. Subject to the three hypotheses, several theorems were derived, for which three model- 
dependent quantities needed to be identified as well: (1) the energy barrier F* separating □ from 
ffl; (2) the set C* of critical configurations; (3) the cardinality N* of the set of protocritical 
configurations, which can be thought of as the "entrance" set of C*. Quantity (1) was identified 
in [15]. In the present paper we identify quantities (2) and (3). 

In [15] the first two hypotheses were verified and the energy value F* of the energy barrier 
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separating □ and ffl is identified. These results were sufRcient to establish the exponential 
probability distribution of the nucleation time divided by its mean, and to determine the mean 
nucleation time up to a multiplicative factor K of order 1 + o(l) as the inverse temperature 

/3 -j> 00. 

Present paper. In the present paper we show that the model satisfies the third hypothesis, 

and we identify the set of critical configurations. We give a geometric characterization of the 
configurations in C* and compute the value of N* . A prototype of the critical configuration was 
already identified in [15], and consists of a configuration of minimal energy with — 1) + 1 
particles of type 2 arranged in a cluster of minimal energy plus a particle of type 2. The difficult 
task is to characterize the full set of critical configurations. This part of the analysis uses the 
specific dynamical features of the model, which are investigated in detail in a neighborhood of the 
saddle configurations. This task is carried over by observing that, in the regime < Ai < t/ < A2 
and in configurations of minimal energy, each particle of type 2 is surrounded by particles of type 
1. This allows us to look at configurations of minimal energy not as clusters of single particles, but 
as clusters of "tiles": particles of type 2 surrounded by particles of type 1 The tiles making up the 
cluster can travel around the cluster faster than particles of type 2 can appear at the boundary 
of A. This motion of tiles along the border gives the dynamics the opportunity to extend the set 
of critical configurations. Different mechanisms are identified that allow tiles to travel around 
a cluster. The energy barrier that must be overcome in order to activate these mechanisms is 
determined, and is compared with the energy barrier the dynamics has to overcome in order 
to let a particle enter A. How rich the set of critical configurations is depends on the relative 
magnitude of these barriers. Consequently, the geometry of the critical configurations is highly 
sensitive to the choice of parameters. 

The problem of computing the value A^"^, i.e., the cardinality of the set of protocritical con- 
figurations, is reduced to counting the number of polyominoes of minimal perimeter belonging to 
certain classes of configurations that depend on the values of the parameters Ai and A2. This 
is a non-trivial problem that is interesting in its own right. With these results we are able to 
derive the sharp asymptotics for the nucleation time and to find the entrance distribution of the 
set of critical configurations. 

Results in this paper are derived by using a foliation of the state space according to the 
number of particles of type 1 in A, plus the fact that configurations in C* must satisfy a "gate 
property", i.e., they must be visited by all optimal paths. These results allow us to compute 
sharp asymptotic values for the expected nucleation time. 

Literature. Similar analyses have been carried out both for conservative and non-conservative 
dynamics. For Ising spins subject to Glauber dynamics in finite volume, a rough asymptotics 
for the nucleation time was derived by Neves and Schoiimann [20] (on 1?) and by Ben Arous 
and Cerf [2] (on 1?). Their results were improved by Bovier and Manzo [7], where the potential- 
theoretic approach to metastability developed by Bovier, Eckhoff, Gayrard and Klein [5[ was 
used to compute a sharp asymptotics for the nucleation time. 

For the model with three-state spins (Blume-Capel model), the transition time and the typical 
trajectories were characterized by Cirillo and Olivieri [8[. For conservative Kawasaki dynamics, 
metastable behavior was studied in den Hollander, Olivieri and Scoppola [13] (on 1?) and in den 
Hollander, Nardi, Olivieri and Scoppola [12] (on 1?). The sharp asymptotics of the nucleation 
time was derived by Bovier, den Hollander and Nardi [6]. Models with an anisotropic interaction 
were considered in Kotecky and Olivieri [17] for Glauber dynamics and in Nardi, Olivieri and 
Scoppola [19] for Kawasaki dynamics. 
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The model studied in the present paper fahs in the class of variations on Ising spins subject 
to Glauber dynamics and lattice gas particles subject to Kawasaki dynamics. These variations 
include staggered magnetic field, ncxt-noarcst-ncighbor interactions, and probabilistic cellular 
automata. In all these models the geometry of the energy landscape is highly complex, and 
needs to be controlled in order to arrive at a complete description of metastability. For an 
overview, see the monograph by Olivieri and Vares [22], Chapter 7. 



Outline. Section 1.1 defines the model. Section 1.2 introduces basic notation and key defini- 
tions, Section 1.3 states the main theorems, while Section 1.4 discusses these theorems. 



1.1 Lattice gas subject to Kawaisaiki dynamics 

Let A C be a large box centered at the origin (later it will be convenient to choose A rhombus- 
shaped). Let I • I denote the Euclidean norm, let 

d-A = {xeA:3y^A:\y-x\ = l}, 

d+A = {x^A:3yeA:\y-x\ = l}, ' 

be the internal, respectively, external boundary of A, and put A~ = A\9~A and A+ = AU 9+A. 
With each site x (E A we associate a variable ri{x) G {0, 1, 2} indicating the absence of a particle 
or the presence of a particle of type 1 or type 2, respectively. A configuration rj = {r]{x) : a; G A} 
is an element oi X = {0,1,2}^. To each configuration r] we associate an energy given by the 
Hamiltonian 

H = -U l^^(^)^(j,)=2} + Ai ^ l{r,(x)=l} + ^2 X! ^Mx)=2}, (1-2) 

(a;,j/)eA*.- xeA xGA 

where A*'~ = {{x,y): x,y e A~, \x — y\ = 1; \x — z\ > 2, \y — z\ > 2 y z e d~ A} is the set of 
non-oriented bonds in A at distance at least 3 from d~A, —U < is the binding energy between 
neighboring particles of different types in A~ , and Ai > and A2 > are the activation energies 
of particles of type 1, respectively, 2 in A. Without loss of generality we will assume that 

Ai < A2. (1.3) 

The Gibbs measure associated with H is 

^^^{^) = ^e-^"^^\ vex, (1.4) 

where /3 e (0, 00) is the inverse temperature and is the normalizing partition sum. 

Kawasaki dynamics is the continuous-time Markov process (r/t)t>o with state space X whose 
transition rates are 

^P(^''^^ = { 0, otherwise, ^^'^^ 

where r] r]' means that r]' can be obtained from ry by one of the following moves: 

• interchanging and 1 or and 2 between two neighboring sites in A 
("hopping of particles in A"), 
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• changing to 1 or to 2 in 5 A 
("creation of particles in d~A"), 

• changing 1 to or 2 to in 9~ A 
("annihilation of particles in 9~A"). 

Note that this dynamics preserves particles in A~ , but allows particles to be created and annihi- 
lated in 9~A. Think of the latter as describing particles entering and exiting A along non-oriented 
bonds between 5+A and d~A (the rates of these moves are associated with the bonds rather than 
with the sites). The pairs (??,?7') with ry -f^ ry' are called communicating configurations, the tran- 
sitions between them are called allowed moves. Note that particles in d~A do not interact: 
the interaction only works well inside A~ (see (1-2)). Also note that the Gibbs measure is the 
reversible equilibrium of the Kawasaki dynamics: 

fJ'0{v)c0iv,v') = fJ'0iv')c0iv',v) V?7,?7'eA'. (1.6) 

The dynamics defined by (1.2) and (1.5) models the behavior in A of a lattice gas in Z^, 
consisting of two types of particles subject to random hopping, hard-core repulsion, and nearest- 
neigbor attraction between different types. We may think of Z^\A as an infinite reservoir that 
keeps the particle densities fixed at pi = e~l^^^, respectively, p2 = e"'^'^^. In the above model 
this reservoir is replaced by an open boundary d~A, where particles are created and annihilated 
at a rate that matches these densities. Thus, the dynamics is a finite-state Markov process, 
ergodic and reversible with respect to the Gibbs measure fj,^ in (1.4). 

Note that there is no binding energy between neighboring particles of the same type (including 
such an interaction would make the model much more complicated). Consequently, our dynamics 
has an "anti-ferromagnetic flavor", and does not reduce to Kawasaki dynamics with one type of 
particle when Ai = A2. Also note that our dynamics does not allow swaps between particles, 
i.e., interchanging 1 and 1, or 2 and 2, or 1 and 2, between two neighboring sites in A. (The first 
two swaps would not effect the dynamics, but the third would; for Kawasaki dynamics with one 
type of particle swaps have no effect.) 



1.2 Bcisic notation and key definitions 

To state our main theorems in Section 1.3, we need some notation. 
Definition 1.1 (a) □ the configuration where A is empty. 

(b) ffl is the set consisting of the two configurations where A is filled with the largest possible 
checkerboard droplet such that all particles of type 2 are surrounded by particles of type 1 . 

(c) w : T] ^ rj' is any (self- avoiding) path of allowed moves from rj G X to rj' G X . 

(d) ^{r],ri') is the communication height between r],r]' G X defined by 

${r),r]') = min maxiJ(^), (1.7) 

and ^{A, B) is its extension to non-empty sets A, B G X defined by 

'P{A,B)= min #(77,7,'). (1-8) 

(e) S{r],r]') is the communication level set between rj and r( defined by 

S{ti,v) = {c e a-: 3w: rj ^ t^' , u B (: maxi?(0 = H{C) = $(ry,7?')} • (1-9) 



5 



A configuration € S{r],r]') is called a saddle for (r?, ??')• 

(f ) is the stability level ofrj&X defined by 

Vr,=<P{v,Irj)-H{n), (1.10) 

where I,, = G : H{^) < H{ri)} is the set of configurations with energy lower than rj. 

(g) -^stab = {r} € X: H{ri) = min^^x H{^)} is the set of stable configurations, i.e., the set of 
configurations with mininal energy. 

(h) A'meta = {t] & X : = ^^iex\Xsta.h ^} of metastable configurations, i.e., the set 
of non-stable configurations with maximal stability level. 

(i) r = for 1] e Xjneta. (notc that r] is constant on X^eta), T* = ^(D, ffl) — H{D) (note 
that H{a) =0). 

Definition 1.2 (a) (77 — >■ ?7')opt is the set of paths realizing the minimax in ^{r],ri'). 

(b) A setW C X is called a gate forr] r( ifW C S{r],T]') andujDW ^ for all to € (?? »7')opt- 

(c) A set W C X is called a minimal gate for rj ^ rj' if it is a gate for rj ^ rj' and for any 
W C W there exists an lo' G {t] r/^opt such that w' n W = 0. 

(d) A priori there may be several (not necessarily disjoint) minimal gates. Their union is denoted 
by Q{ri, ij') and is called the essential gate for [rj — > 77')opt- The configurations in S{r], ri')\Q{rj, r]') 
are called dead-ends. 

(e) Let S{w) = {argmax^g;^ A saddle ( S S{i], rj') is called unessential if, for all co e {ij ^ 
il')opt such that a; 3 C the following holds: S{lo)\{(} ^ and there exists an w' e (r? — )■ ?7')opt 
such that S{uj') C 5'(w)\{C}. 

(f) A saddle ( £ S{r], ri') is called essential if it is not unessential, i.e., if either of the following 
occurs: 

(fl) There exists an co G {ri ^ v')opt such that S{oj) = {Q}. 

(£2) There exists an lv G (?? ^ V)opt such that S{lo) 3 {Q and S{lo') ^ S'(a;)\{C} for all 
w' e (7? 7?')opt- 

Lemma 1.3 [Manzo, Nardi, Olivicri and Scoppola [18], Theorem 5.1] 
A saddle C, S S{r],r]') is essential if and only if C ^ G(j]tV')- 

In [14] we are interested in the transition of the Kawasaki dynamics from □ to ffl in the hmit 
as ^ — ^ 00. This transition, which is viewed as a crossover from a "gas phase" to a "hquid phase", 
is triggered by the appearance of a critical droplet somewhere in A. The critical droplets form a 
subset C* of the essential gate G{0, ffl), and all have energy F* (because H(0) = 0). 

In [14] we showed that the first entrance distribution on the set of critical droplets is uniform, 
computed the expected transition time up to and including a multiplicative factor of order one, 
and proved that the nucleation time divided by its expectation is exponentially distributed, all 
in the limit as /3 ^ cjo. These results, which are typical for metastable behavior, were proved 
under three hypotheses: 

(HI) -Ystab = ffl- 

(H2) There exists a < such that Vr, < V* for all j? e X\{a, ffl}. 
(H3) See (H3-a,b,c) and Fig. 1 below. 

The third hypothesis consists of three parts characterizing the entrance set of Q{0, ffl) and the 
exit set of G{0, ffl). To formulate these parts some further definitions are needed. 
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Definition 1.4 (a) C^^ is the minimal set of configurations in Q{D,S) such that all paths in 
(□ -> ffl)opt enter 0(0, ffl) through C^^. 

(b) V is the set of configurations visited by these paths just prior to their first entrance of 5(0, ffl). 

(H3-a) Every fj d V consists of a single droplet somewhere in . This single droplet fits inside 
an X square somewhere in for some L* G N large enough that is independent of 
fj and A. Every rj G C^^ consists of a single droplet r) € V and one additional free particle 
of type 2 somewhere in d^A. 

Definition 1.5 (a) C*^^ is the set of configurations obtained from C^^ by moving the free particle 
of type 2 along a path of empty sites in A and attaching it to the single droplet (i.e., creating at 
least one additional active bond). This set decomposes as C*^^ — Ufi^-pC*^^{ri) . 
(b) C* is the set of configurations obtained from C^^ by moving the free particle of type 2 along 
a path of empty sites in A. This set decomposes as C* — Ufj^-pC*{rj). 

Note that F* = H{C*) = H{V) + A2, and that C* consists of precisely those configurations 
"interpolating" between V and C*tj: a free particle of type 2 enters i9~A and moves to the single 
droplet where it attaches itself via an active bond, i.e., a bond between particles of type 1 and 
2. Think of V as the set of configurations where the dynamics is "almost over the hill", of C* as 
the set of configurations where the dynamics is "on top of the hill", and of the free particle as 
"achieving the crossover" when it attaches itself "properly" to the single droplet (the meaning of 
the word "properly" will become clear in Section 5; see also [14], Section 2.4). The sets V and 
C* are referred to as the protocritical droplets, respectively, the critical droplets. 

(H3-b) All transitions from C* that either add a particle in A or increase the number of droplets 
(by breaking an active bond) lead to energy > F*. 

(H3-c) All oj e (C|^j E)opt pass through C*^^. For every fj e V there exists a ^ G Qtt('?) such 



Figure 1: A qualitative representation of a configuration in C^^. If the free particle of type 2 
reaches the site marked as then the dynamics has entered the "basin of attraction" of ffl. 

Remark: Hypothesis (H3-a) and Definition 1.5 are slightly different from how they appear in 
[14] and ]15]. This is done to make their analogues in ]14] and [15] more precise, and to allow for 
a more precise proof of Lemma 1.18 and Lemma 2.2 in [14], which we repeat in Appendix B.l. 

As shown in [14], (H1-H3) are needed to derive the metastability theorems in [14] with the 
help of the potential-theoretic approach to metastability outlined in Bovier ]3]. In ]15] we proved 



that $(C,ffl) < F*. 



L* 




L* 
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(H1-H2) and identified the energy F* of critical droplets. In the present paper we prove (H3), 
identify the set C* of critical droplets, and compute the cardinality N* of the set V of protocritical 
droplets modulo shifts, thereby completing our analysis. 

Hypotheses (H1-H2) imply that (A'meta, -^stab) = and that the highest energy barrier 

between a configuration and the set of configurations with lower energy is the one separating 
□ and ffl, i.e., (Diffl) is the unique metastable pair. Hypothesis (H3) is needed to find the 
asymptotics of the prefactor of the expected transition time in the limit as A ^ and will 
be proved in Theorem 1.7 below. The main theorems in [14] involve three model-dependent 
quantities: the energy, the shape and the number of critical droplets. The first (F*) was identified 
in [15], the second (C*) and the third {N*) will be identified in Theorems 1.8-1.10 below. 



1.3 Main theorems 

In [14] it was shown that < A1+A2 < 4t/ is the metastable region, i.e., the region of parameters 
for which □ is a local minimum but not a global minimum of H. Moreover, it was argued that 

within this region the subregion where Ai, A2 < [/ is of little interest because the critical droplet 
consists of two free particles, one of type 1 and one of type 2. Therefore the proper metastable 
region is 

< Ai < A2, Ai + A2 < 4U, A2 > U, (1.11) 

as indicated in Fig. 2. 



Ai + Aj = 4(7 




Figure 2: Proper metastable region. 



In this present paper, as in [15], the analysis will be carried out for the subregion of the 
proper metastable region defined by 

Ai <U, A2 - Ai > 2U, Ai + A2 < 4U, (1.12) 

as indicated in Fig. 3. {Note: The second and third restriction imply the first restriction. 
Nevertheless, we write all three because each plays an important role in the sequel.) 

Hypothesis (H3) involves additional characterizations of the sets V and C*. It turns out that 
these sets vary over the region defined in (1.12). The subregion where A2 < 4C/ — 2Ai is trivial: 
the configurations in V consist of a single droplet, with one particle of type 2 surrounded by four 
particles of type 1, located anywhere in A~. For this case, F* = 4Ai + 2A2 — iU and A''* = 1. 
We will split the subregion where A2 > 4:U — 2Ai into four further subregions (see Fig. 4). 

For three of these subregions we will indentify V, C* and C^^, prove (H3), and compute N*, 
namely. 
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Figure 3: Subregion of the proper metastable region given by (1.12). 




U Ai 



Figure 4: Subrcgions of the parameter space. In the black region: I* < 3. The regions RA, RB 
and RC are, respectively, light gray, dark grey and dashed. 
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RA: Ai < 3U; 

RB: 3C/ < A2 < 2;7 + 2Ai; 

RC: A2 > 3C/ + Ai. 

The fourth subregion is more subtle and is not analyzed in detail (see Section 1.4 for connnents). 
All subregions are open sets. This is done to avoid parity problems. We also require that 



Ai/e^N with £ = 4C/ - Ai 



and put 



Ai 



eN\{i}. 



(1.13) 



(1.14) 



To state our main theorem we need the following definitions. A 2-tile is a particle of type 2 
surrounded by four particles of type 1. Dual coordinates map the support of a 2-tile to a unit 
square. A monotone polyomino is a polyomino whose perimeter has the same length as that 
of its circumscribing rectangle. Given a set of configurations V, we write V^^"^ to denote the 
configurations obtained from V by adding a particle of type 2 to a site in d~A. (For precise 
definitions see Sections 2.1-2.2.) 



Definition 1.6 (a) Va is the set of2-tiled configurations with £*{£* — 1) + 1 particles of type 2 
whose dual tile support is a rectangle of side lengths £*,£*— 1 plus a protuberance on one of the 
four side of the rectangle (see Fig. 14 

(b) Vb is the set of 2-tiled configurations with £*{£* — 1) + 1 particles of type 2 whose dual 
tile support is a monotone polyomino and whose circumscribing rectangle has side lengths either 
£*,£* or£* + 1,£* - 1 (see Fig. 2Q). 

(c) Vc is the set of 2-tiled configurations with £*{£* — 1) + 1 particles of type 2 whose dual 
tile support is a monotone polyomino and whose circumscribing rectangle has perimeter (see 
Fig. 27). 

Note that VaCVb CVc. 

Theorem 1.7 Hypothesis (H3) is satisfied in each of the subregions RA — RC. 
Theorem 1.8 In subregion RA, V = Va, C^^ = V^'^, and N* = 8£* - 4. 
Theorem 1.9 In subregion RB, V = Vb, C^d = ^* = + ri*^i-i]. 



Theorem 1.10 In subregionRC, V = Vc, C^j = ^c*^ andN* = 8[qe*-i+Y^[l^^ r^._c2_i]. 

Here, (r^) and {qk) are the coefficients of two generating functions defined in Appendix A, which 
count polyominoes with fixed volume and minimal perimeter. The claims in Theorems 1.8-1.10 
are valid for > 4 only. For i* = 2, 3, see Section 4.3. 



1.4 Discussion 

1. In (1.2) we take an annulus of width 3 without interaction instead of an annulus of width 1 
as in [14] and [15]. This allows us to prove that the model satisfies (H3), without having to deal 
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with complications that arise when droplets are too close to the boundary of A. The theorems 
in [14] and [15] remain valid. 

2. Theorems 1.7 and 1.8-1.10, together with the theorems presented in [14] and [15], complete 
our analysis for part of the subregion given by (1.12). Our results do not carry over to other 
values of the parameters, for a variety of reasons explained in [14], Section 1.5. In particular, for 
Ai > [/ the critical droplets are square-shaped rather than rhombus-shaped. Moreover, (H2) is 
expected to be much harder to prove for A2 — Ai < 2U. 

3. Theorems 1.8-1.10 show that, even within the subregion given by (1.12), the model-dependent 

quantities C* and A^*, which play a central role in the metastability theorems in [14], are highly 
sensitive to the choice of parameters. This is typical for metastable behavior in multi-type 
particle systems, as explained in [14], Section 1.5. 

4. The arguments used in the proof of Theorems 1.7 and 1.8-1.10 are geometric. Along any 
optimal path from □ to ffl, as the energy gets closer to F* the motion of the particles becomes 
more resticted. By analyzing this restriction in detail we are able to identiiy the shape of the 
critical droplets. 

5. The fourth subregion is more subtle. The protocritical set V is somewhere between T>b 
and Vc, and we expect V = Vc for small Ai and Vb ^ V C. Vc for large Ai. The proof of 
Theorems 1.8-1.10 in Section 5 will make it clear where the difficulties come from. 

Outline: In the remainder of this paper we provide further notation and definitions (Section 2), 
state and prove a number of preparatory lemmas (Section 3), describe the motion of "tiles" along 
the boundary of a droplet (Section 4), and give the proof of Theorem 1.7 and 1.8-1.10 (Section 5). 
In Appendix A we recall some standard facts about polyominoes with minimal perimeter. 



2 Coordinates and definitions 

Section 2.1 introduces two coordinate systems that are used to describe the particle configura- 
tions: standard and dual. Section 2.2 lists the main geometric definitions that are needed in the 
rest of the paper. 



2.1 Coordinates 

1. A site i G A is identified by its standard coordinates x{i) = {xi{i),X2{i)), and is called odd 
when xi{i) + X2{i) is odd and even when xi{i) + X2{i) is even. The standard coordinates of a 
particle p in A are denoted by x{p) = (.Ti(p), X2{p))- The parity of a particle p in a configuration 
77 is defined as X\{p) + X2{p) + r]{x{p)) modulo 2, and p is said to be odd when the parity is 1 
and even when the parity is 0. 

2. A site i e A is also identified by its dual coordinates 

^^^^ ^ xAi)_X2(^^ ^^^^ ^ x,{^)+X2{^) _ ^^.1) 

Two sites i and j are said to be adjacent, written i j, when \x\{i) — x\{j)\ + \x2{i) — X2{j)\ = 1 
or, equivalently, \ui{i) - ui(j)| = \u2{i) - W2(i)| = \ (see Fig. 5). 
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3. For convenience, we take A to be the x dual square with bottom-left corner at site 

with dual coordinates (--^, --^) for some i e N with L > 2i* (to allow for H{m) < H{n)). 
Particles interact only in A~, which is an (L + i) x (L + ^) dual square. This dual square, 
a rhombus in standard coordinates, is convenient because the local minima of H are rhombus- 
shaped as well (for more details see [15]). 




(a) (b) 



Figure 5: A configuration represented in: (a) standard coordinates; (b) dual coordinates. Light- 
shaded squares are particles of type 1, dark-shaded squares are particles of type 2. In dual 
coordinates, particles of type 2 are represented by larger squares than particles of type 1 to 
exhibit the "tiled structure" of the configuration. 



2.2 Definitions 

1. A site i S A is said to be lattice- connecting in the configuration rj if there exists a lattice 
path A from i to d~A such that ri{j) = for all j G X with j ^ i. We say that a particle p is 
lattice-connecting if x{p) is a lattice-connecting site. 

2. Two particles in 77 at sites i and j arc called connected if i ^ j and ri{i)r]{j) = 2. If two 
particles pi and p2 are connected, then we say that there is an active bond b between them. 
The bond b is said to be incident to pi and P2- A particle p is said to be saturated if it is 
connected to four other particles, i.e., there are four active bonds incident to p. The support 
of the configuration rj, i.e., the union of the unit squares centered at the occupied sites of 77, 
is denoted by supp(?7). For a configtiration jj, ni{r]) and n2{ri) denote the number of particles 
of type 1 and 2 in 77, and -6(77) denotes the number of active bonds. The energy of 77 equals 
H{r]) = Amiirt) + A2n2{v) - UB{r]). 

3. Let G{r]) be the graph associated with 77, i.e., G{r]) = {V{r]), E{r])), where 1^(77) is the set of 
sites i G A such that ri{i) ^ 0, and E{ri) is the set of pairs {i,j}, i,j G V{r]), such that the 
particles at sites i and j are connected. A configuration r/' is called a subconfiguration of rj, 
written 77' -< 77, if 77' (i) = 77(i) for alH € A such that 77' (i) > 0. A subconfiguration c ^ 77 is called 
a cluster if the graph G{c) is a maximal connected component of G{i]). The set of non-saturated 
particles in c is called the boundary of c, and is denoted by dc. Clearly, all particles in the same 
cluster have the same parity. Therefore the concept of parity extends from particles to clusters. 

4. For a site 7 € A, the tile centered at i, denoted by t(7), is the set of five sites consisting of i 
and the four sites adjacent to i. If i is an even site, then the tile is said to be even, otherwise the 
tile is said to be odd. The five sites of a tile are labeled a, b, c, d, e as in Fig. 6. The sites labeled 
a, b, c, d are called junction sites. If a particle p sits at site i, then t{i) is alternatively denoted 
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by t{p) and is called the tile associated with p. In standard coordinates, a tile is a square of size 
\/2. In dual coordinates, it is a unit square. 



5. A tile whose central site is occupied by a particle of type 2 and whose junction sites are 
occupied by particles of type 1 is called a 2-tile (see Fig. 6). Two 2-tiles are said to be adjacent 
if their particles of type 2 have dual distance 1. A horizontal (vertical) 12-bar is a maximal 
sequence of adjacent 2-tiles all having the same horizontal (vertical) coordinate. If the sequence 
has length 1, then the 12-bar is called a 2-tiled protuberance. A cluster containing at least one 
particle of type 2 such that all particles of type 2 are saturated is said to be 2-tiled. A 2-tiled 
configuration is a configuration consisting of 2-tiled clusters only. A hanging protuberance (or 
hanging 2-tile) is a 2-tile where three particles of type 1 are adjacent to the particle of type 2 
of the 2-tile only (see Fig. 17(b)). 

Remark 2.1^4 configuration consisting of a dual 2-tiled square of side length i* belongs to X^. 























(a) 



(b) 



H 





(c) 



(d) 



Figure 6: Tiles: (a) standard representation of the labels of a tile; (b) standard representation 
of a 2-tile; (c) dual representation of the labels of a tile; (d) dual representation of a 2-tile. 



6. The tile support of a configuration rj is defined as 

M= U ^(P)' (2.2) 

p£ro2(77) 

where W2 (rj) is the set of particles of type 2 in 77. Obviously, [rj] is the union of the tile supports 
of the clusters making up rj. For a standard cluster c the dual perimeter, denoted by P(c), is the 
length of the Euclidean boundary of its tile support [c] (which includes an inner boundary when 
c contains holes). The dual perimeter P(r?) of a 2-tiled configuration 7/ is the sum of the dual 
perimeters of the clusters making up r/. 

7. Denote by V*,n2 the set of configurations such that in (A~)~ the number of particles of type 

2 is n2. Denote by Vf"^^ the subset of V*,n2 where the number of active bonds is 4n2 and there 
are no non-interacting particles of type 1, i.e., the set of 2-tiled configurations with 712 particles 
of type 2. A configuration 77 is called standard if 77 G V^"^^ ^^'^ its tile support is a standard 
polyomino in dual coordinates (see Definition 2.2 below). A configuration 77 with 712(77) particles 
of type 2 is called quasi-standard if it can be obtained from a standard configuration with 712 (77) 
particles of type 2 by removing some (possibly none) of the particles of type 1 with only one 
active bond, i.e., corner particles of type 1. Denote by V*,„ the set of configurations of minimal 
energy in V*,„. 

8. The state space X can be partitioned into manifolds: 

^= [j V.,„2- (2.3) 

n2=0 
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Two manifolds V*,„ and V*,„' are called adjacent if |n— n'| = 1. Note that transitions between two 
manifolds are possible only when they are adjacent and are obtained either by adding a particle 
of type 2 to 9~A (V+.„ V*.„+i) or removing a particle of type 2 from d~A (V+,„ — > V*,„_i). 
Note further that □ S V*,o and ffl S V*,l2. Therefore, to realize the transition □ — >• ffl, the 
dynamics must visit all manifolds V*,„ with n = 1, . . . , L^. Abbreviate V*,<to = [Xi=o ^*,n and 

9. For y C X, fj € y and x G A\supp[r7], we write r] = {fj,x) to denote the configuration that 
is obtained from fj by adding a particle of type 2 at site x. We write y^'^'^ to denote the set of 
configurations obtained from a configuration in y by adding a particle of type 2 in d~A, i.e., 
•ybd2 _ \J^^y\J^^Q_^{fj,x). For w. □ ^ ffl, let <Ty{u)) be the configuration in y that is first 
visited by u. Define 

fy= U M^), 

optimal 

called the entrance, respectively, the optimal extrance of y. With this notation we have C^^ = 

10. For A,BCX, define 

g{A,B) = {7jeB: 3 ( e A '.md u; : ( ^ tj: n^iv) < n2{0 < n^iC), H{0 <T*y ^ € u;}, 
giA,B) = {r] € B: 3C € A and u: C ^ r]-. n2{v) < ^2(0 < MO, H{^) < r*V^ €uj}. ^ ' ^ 



(a) (b) (c) 

Figure 7: Corners of polyominoes: (a) one convex corner; (b) one concave corner; (c) two concave 
corners. Shaded mean occupied by a unit square. 

11. A unit hole is an empty site such that all four of its neighbors are occupied by particles of 
the same type (either all of type 1 or all of type 2). An empty site with three neighboring sites 
occupied by a particle of type 1 is called a good dual corner. In the dual representation a good 
dual corner is a concave corner (see Fig. 7). The surface of 77 € A* is defined as 

F(r?) = {xeA:3y^x: r]{y) = 1}. (2.6) 

For ri€ X, let 

r{v) = + - m] = 2P(,7) + 4[C(ry) - Q(r;)], (2.7) 

where C{r]) is the number of clusters in t], P{r]) the total length of the perimeter of these clusters, 
Q{r]) the number of holes, tp{r]) the number of convex corners, and (p{f]) is the number of concave 
corners. Note that T{ri) = '}2c^^T{c), where the sum runs over the clusters in rj. 

We also need the following definition: 
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Definition 2.2 [Alonso and Cerf [1].] A polyomino (= a union of unit squares) is called mono- 
tone if its perimeter is equal to the perimeter of its circumscribing rectangle. A polyomino is 

called standard if its support is a quasi-square (i.e., a rectangle whose side lengths differ by at 
most one), with possibly a bar attached to one of its longest sides. 

3 Preparatory lemmas 

In this section we collect a number of preparatory lemmas that are valid throughout the subregion 
given by (1.12). These lemmas will be needed in Section 5 to prove Theorems 1.7 and 1.8-1.10. 
In Section 3.1 we characterize [*V*^^*(^*_i)+i (Lemmas 3.2-3.3 below), in Section 3.2 we charac- 
terize (;({□}, V*_£.(£._i)_|_i) and (;({□}, V^^^*(^*_i)_|_i) (Lemma 3.3 below), and in Section 3.3 we 
characterize Q{D,S) (Lemmas 3.5-3.6 below). 

An elementary observation is the following: 

Lemma 3.1 Ifrj £ V*,n2 with i* {i* — 1) + 1 <n2< {i*)^, thenr] is 2-tiled and its dual perimeter 
is equal to 4£*. 

Proof. Immediate from Lemmas 2.2-2.3 and 4.1 in [15], and also from Corollary 2.5 in [1]. ■ 
3.1 Chciracterization of f^V^t,^*(£*_i)+i 

Lemma 3.2 Let p be a2 tiled configuration with £*{£* — 1) particles of type 2 a,nd with dual tile 
support equal to a rectangle of side lengths (.* , {£* — 1) (i.e., p is a standard configuration). 

(1) If T}^ p is a 2-tiled configuration with £*{£* — 1) particles of type 2, then H{r]) > F* — A2. 

(2) If Tj ^ p is a configuration with £*{£* — 1) particles of type 2 such that [r}] ^ [p], then 
H{ri) >r*- A2. 

(3) Ifrj^pisa configuration with £*{£* — 1) particles of type 2 obtained from p by removing at 
least one of the "non-corner" particles of type 1 in p (note that r] and p have the same dual tile 
support), then H{r]) > F* — A2. 

Proof. (1) Since 77 is a 2-tiled configuration, it follows from Lemma 2.3 in [15] that II{rf) — H{p) = 
\\T{r}) — T(p)]Ai, because the energy difference between the two configurations only depends 
on the difference in the number of particles of type 1. From Lemma 2.2 in [15] it follows that 
T{ri) > T{p). From the definition of T in (2.7) and Eq. (2.4) in [15] we have that, for any 2-tiled 
r], T{ri) = 4k for some k € N. Hence T{r]) - T{p) > 4, and so H{r]) - H{p) > Ai. The claim 
now follows by observing that II{p) = F* + £ — Ai — A2 and e > 0. 

(2) First consider the case where 77 consists of a single cluster. Then there exists a configuration 
rj' , obtained from rj by saturating all particles of type 2, such that II{r]') < H{ri) with equality 
if and only if r}' = rj. Clearly, [r]] = [r]']. By part (1), we have II{ri) > H{ri') > F* — A2. If 77 
consists of clusters ci, . . . ,Cm with m G N\{1}, then observe that i?(ry) — Y^T=i '^i- ?y"2(ci) 
denote any standard configuration with n2(ci) particles of type 2. By Lemmas 3.1 and 4.2 in 
[15], we have II{ri) — Yli=i-^i^'i) — 5^iLi -^('7"^'"^''')- Since p is a standard configuration, it 
follows from Lemma 2.2 in [15] that X]i=i '^('7"^'''^''') > T{p)j ^^'^ ^o, as in the proof of part 
(1)' '^('?"^*''^''') ~ Tip) > 4. Using (3.6) in [15] for the energy of a standard configuration, 

we obtain that Yli=i H{rf^(''*"i) - H{p) = ^[ELi T(?7"=(^*)) - T(p)]Ai, from which we get the 
claim. 
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(3) Let m G N denote the number of non-corner particles of type 1 removed from p to obtain rj. 
Then H{r]) > H{p) + m{2U — Ai) > H{p) + 2U — Ai (because each of the non-corner particles of 
type 1 in p has at least 2 active bonds). Substituting the value of H{p) into the latter expression, 
we obtain ff (r?) > T* - Ai - A2 + £ + 2?/ - Ai . The claim follows by observing that Ai < J7. ■ 

Lemma 3.3 (1) All paths in (□ ffl)opt enter the set V*,£*(£*_i)_,_i via a configuration {fj,x) 

with r) G Vi,^£.(£._i) a quasi- standard configuration and x € d~ K. 

(2) All paths in (□ — ffl)opt enter the set V^^£*(£*_i)_j_2 via a configuration {r),x) with f] G 
V^^^*(^*_i)+i such that ^{D,fj) < T*, i.e., fj G V*,£*(^*-i)+i), andx G d~A. Consequently, 

^({□}, Vj,_^.(£*_i)+i)^'^^ is a gate for the transition □ — )• ffl. 

Proof. (1) This is immediate from Lemma 3.2. 

(2) By Theorem 1.5 in [15] (which identifies F*) and Lemmas 3.1-3.2 in [15] (which determine 
the energy of configurations in V*.„ for all n), if 77 € Vj,^^.(£*_i)_|_i, then H{r)) = F* — A2. We 
argue by contradiction. Suppose that cj G (□ — > ffl)opt enters Vi,^£*(£*_i)+2 via a configuration 

C = {C,x) with C e V^,^*(^*_i)+i\K,^*(£*-i)+i and X S 9" A. Then H{C) = H{C) + A2 > F*, 
because H{Q > F* — A2. Hence u is not optimal. ■ 

3.2 Characterization of 5f({n}, and ^({D}, V*/*(<?*_i)+i) 

Definition 3.4 (a) For n G N, let Sn be the set of standard configurations with n particles of 
type 2. 

(b) Let uj: □ ^ EE = (□, . . . , ^, 77, C, • • • , fflj- Write Pi^{ri) to denote the part ofu from to ^ and 
Suiiv) to denote the part of co from to W. Any configuration in Pu{r]) is called a predecessor 
of r] in UJ, while any configuration in Saj{ri) is called a successor ofrj in uj. The configurations ^ 
and C are called the immediate predecessor, respectively, the immediate successor ofrj in uj. 

Lemma 3.5 (1) For every ( G 9{{^}j'^*,e*(e*-i)+i) there is a standard configuration fj G 

St;*^i>*-i)+i such that ( G g({?7}, V*,£*(£*_i)+i). 

(2) For every ( G (;({□}, Vi,.^*(^*_i)+i) there is a standard configuration fj G «S£*(^._i)_|_i such 
that C e Consequently, 

»)65f*(^*_i)+i 

Proof (1) Pick C e .9({n}, V*.^*(f._i)+i). Let u: □ ^ C be such that max^^^ i7(^) < F*. Let 
r] be the configuration visited by oj when it enters the set V*,^*(£*_i)+i for the last time before 
visiting (. Write u) as wi +0J2, where uji is the part of w from □ to 77 and UI2 is the part of w from 77 
to By Lemma 3.2, we have 77 = (jj,x), where 77 is a quasi-standard configuration in Vi, £.(£._i) 
and X G d~A, otherwise -^(77) > F*. We will show that there is a standard configuration 
?7 G V*,^*(£*_i)_|_i and a path W3: r] ^ fj such that H{^) < H{ri) and 7i2(^) = £*{£* — 1) + 1 for 
all ^ G UJ3. 

Let rj be the standard configuration in Vy,.(^*(^*_i) with the same tile support as rj. This 
configuration exists because every quasi-standard configuration whose support lies in has no 
particle of type 2 in d~A~. (The latter is due to the fact that, in a quasi-standard configuration 
with — 1) particles of type 2, each site that is occupied by a particle of type 2 has at least 
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three neighboring sites occupied by a particle of type 1, and ah sites in d~A~ have at most 
two adjacent sites in A~.) Let 77 the standard configuration in V*^^* (£._!)+! obtained from 77 by 
adding a protuberance, with the particle of type 2 in this protuberance located at a site y* on 
one of the longest sides of the rectangular cluster of 77. This is always possible because at least 
one of the longest sides of [f]] is far away from d~A. 

Consider the path S^^iv)- Since 77 is the configuration visited by to when the set Vi,_£*(£*_i)+i 
is entered for the last time before visiting all configurations in Saj2 {v) have at least — 1) + 1 
particles of type 2. In particular, the particle of type 2 in a; cannot leave A. We refer to this 
particle as the "floating particle". Observe that ^^(77) > F* + e — Ai, with equality if and only 
if 77 is standard. This implies that only moves of the floating particle are allowed until it enters 
A~ (particles in d~A cannot have active bonds). Furthermore, since L > 2£* (and hence the 
sides of fj are smaller than the sides of A) , it follows that all sites y G A~ such that y ^ supp(?7) 
are lattice-connecting. In particular, there exists a lattice path A = a;o, a;i, . . . , Xm in A for some 
7n e N with xq = X and Xm = y*- 

Let W3 be the path from rj to fj obtained by flrst letting the floating particle move along the 

lattice path A until it reaches site y* and then saturating all the particles of type 2 in {fi,y*). 
Note that H{fj,y*) < H(r]) — U and that all conflgurations in 73 have at least as many active 
bonds as 77. Therefore H{^) < H{ri) for all ^ e W3. 

Let 0)3 denote the path from fj to rj obtained by inverting ^3. Then, by construction, Q = 
W3 + W2 is a path from 77 to C such that H{^) < T* and 7x2(0 ^ ^*(-^* — 1) + 1 for aU ^ e d). 

(2) Same as part (1). ■ 



3.3 Characterization of g{D, ffl) 

Lemma 3.6 (1) AU saddles m Vj,^<£*(£*_i) are unessential. 

(2) Let C e V^^<t*{e*-i)+i be suchl.hat H{C) = T* . Let 0(0 = {tJ G (□ ^ ffl)opt : uj ^ (,} he 
the set of optimal paths visiting Q. If all paths in 0{C,) visit 5({n}, V*,£*(£*-i)+i) after visiting 
C, then C is unessential. 

Proof. (1) Let Q S V*,<t*(t*-i) be a configuration such that H{Q) = T* . By Lemma 3.3(2), we 
have (recall Definition 1.2(e)) S'(a;)\{C} 7^ for all w G 0{(), i.e., all paths in 0(C) visit at least 

one other saddle configuration. Pick w €E (□ ^ ffl)opt- By Lemma 3.3(1), uj (last) enters the set 
V^^f*(f*_i)+i via a configuration {pi, x) with pi a quasi-standard configuration obtained from the 
standard configuration pQ G V*^£*(f*_i) by removing i corner particles of type 1, and x G d~A. It 
is clear that the configuration visited by uj just before {pi, x) is pi. Write uj = uji +uj2, where uji is 
a path from □ to pi and W2 is a path from pi to ffl. Obviously, ( G uii. Moreover, w e (□ ffl)opt 
implies that H{pi,x) < T* and, consequently, H{pi) < F* — A2. Furthermore, H{pj) < H{pi) 
for j < i. Let lu^: □ — > po be a path from □ to po such that -ff (0 < T* for all ^ G uj^ (e.g. follow 
the construction of the "reference path" in [15]), and let lo^: po ^ pi be the path obtained by, 
iteratively, detaching and moving out of A one corner particle of type 1 until configuration pi is 
reached. It is easy to see that max^^^^ < F*. Consider the path uj = uj^ + UJ4 + u>2- By 

construction, w G (□ — > E8)opt and 5(a)) C 5'(w)\{C}- Finally, observe that the same argument 
holds for any uj G 0{(). 

(2) Since g{{fj},Vi,^e*(^e*-i)+i) C V*,£*(£*_i)+i, all optimal paths from g{{f]},%j*(^e*_i)+i) to ffl 
must visit the set ^Vi^/*(^e*-i)+2- Therefore S'(w)\{C} 7^ for all uj € C(C): since all configu- 
rations in r^Vy,_^*(^*_i)+2 have energy F*. Pick uj G C(C), and let 77 be the first configuration 
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in <?({»?}, Vjf_£*(^*_i)+i) visited by w after visiting (. Let uii be the part of ui from D to r] and 
u)2 the part of u from 77 to ffl. Since r] € ^({j?}, V*.£*(£*-i)+i), there is a path ^3 from □ to 77 
such that H{a) < F* for all a E ui^. Let di = uj^ + ^2- By construction, a) G (□ — >■ ffl)opt and 
S{w) C S'(w)\{C}. Finally, observe that the same argument holds for all w S C(C)- ' 

Lemma 3_.7 V*/*(r-i)+i) = g-C^, then C^^ C fif(n, anrf 

P C S'(n,V*,£.(£._i)+i). 

Proo/. By Lemma 3.3(2), if 5(0, V^,,*(,*_i)+i) =_(/(□, then <;(□, 

is a gate for the transition □ — >■ ffl. Since (?(□, V*,£*(£*_i)+i)'^'^*^ is a gate, there exists a W C 
(/(□, V*.£*(£*_i)+i)'"^^ that is a minimal gate for the transition □ ^ ffl. Let cj e (□ ^ ffl)opt- 
Since W C g{D,a) and w n W ^ 0, it follows that w_n ^?(D,ffl) n c/(D, ^ 0. 

Combining Lemmas 3.3 and 3.6, it follows that if ^(D, Vjf^^*(^*_i)+i) = g'(n, Vjf^^*(^*_i)+i), then 
all saddles in V*^<£*(£*_i)+i are unessential and, by Lemma 1.3, do not belong to ^(□,ffl). 
Therefore the first configuration in g{0, ffl) visited by ui is an element of g{D, V^,,£*(^*_i)_|_i)'"^^. 
Since the choice of a; e (□ ^ ffl)opt is arbitrary, we conclude that Cj!^^ c g{0, V^,_£*(£*_i)_|_i)'^^^. 

It remains to show that V C ^(n, V*,^*(^._i)+i). The proof is by contradiction. Pick 
fj e P\5'(n,K,£*(<>._i)+i). Since C^^ C ,g(n, K,i^*(i^*_l)+l)'''^^ there is a configuration r] G 
^(n, V,t_£*(£*_i)+i)''''^ obtained in a single step from r}. Clearly, i) G Vy,_>£*(£*_i)_|_i, since the 
number of particle of type 2 in A changes at most by one at each step. Since, by Lemma 3.3 
and the hypothesis g{D,V^ f*(f*_i^_^_i) = .g(n, Vj,.£*(£*_i)_|_i), all paths in (□ ffl)opt enter 
Vi,^e*{e*-i)+2 by adding a particle of type 2 in d~A to a configuration in g(\3,Vi,^i*(^e*-i)+i), it 
follows that fj e V*,>£*(£*_i)+2 by assumption. ■ 



4 Motion of 2-tiles 

In Section 4.1 we study the motion of 2-tiles. In Section 4.2 we derive some restrictions on 
the transitions between configurations with different tile support. In Section 4.3 we identify the 
critical droplets for small values of £*, namely, ^* = 2, 3. 

4.1 Motion of dimers of 2— tiles 

Definition 4.1 (a) Two configurations rj and-q' (with the same tile support) are called equivalent 
if there is a path uj: r] ^ rj' (possibly of length zero) such that all configuration in uj have the 
same tile support and <P{r],r]') < F*. In words, two configurations are equivalent if it is possible 
to go from one to the other via a sequence of moves of particles of type 1 without reaching energy 
level T*. 

(b) A heavy-step is a sequence of moves realizing the transition between two configurations ri and 
rf with different tile support. Note that a heavy-step is completed by moving, removing or adding 
a particle of type 2 in A. 

Let rj G V*,„2 and fj G V*,„2- By Lemma 4.1 and the proof of Lemma 3.1 in [15], both B{fj) 
and ni {fj) are constant in V*,„2 • 

Definition 4.2 (a) A configuration rj G Vi,^n2 "^^ to have m broken bonds if B(ri) = B(fj) — m 
for all fj G V*,n2 • The number of broken bonds in configuration r] is denoted by B~ {rf) . 
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(b) A configuration rj € V*,„2 have n extra particles of type 1 if n\{r]) = ni{fj) — n for 
all f] e V*,„2 . The number of extra particles of type 1 in configuration r] is denoted by nf (rj) . 

(c) B(p, rj) denotes the number of active bonds adjacent to particle p in configuration rj. 

(d) A dimer consists of two adjacent particles of different type such that the particle of type 1 is 
lattice- connecting and has only one active bond (i.e., is a comer particle of type 1). The particle 
of type 2 belonging to a dimer is called a corner particle of type 2. 

(e) A particle of type 2 in a 2-tiled configuration rj is called external if it can be moved without 
moving any other particle of type 2 inr] (see Fig. lA). 

In this section we will exhibit two methods to move a dimer in a configuration to a good 
dual corner of the cluster it belongs to (see Fig. 8). The configuration r]' that is obtained in this 
way satisfies H{ri') < H{ri). In particular, we will exhibit two different choices for a path lu from 
T] as in Fig. 8(a) to rj' as in Fig. 8(b), and we will determine max^g^^ In what follows we 

write AH{co) = H{fi) — H{ri), where f] is the configuration that is reached after the last step in 

CO. 




(a) (b) 



Figure 8: Motion of a 2-tile. 

Lemma 4.3 A 2-tile can be moved within energy barrier 3U and U + 4Ai. 

Proof. We will give examples that are paradigmatic for the general case. Let p, q denote, 
respectively, the particle of type 2 and of type 1 of the dimer that we want to move. 

1. The first method is achieved within energy barrier 3U energy (i.e., max^g^^ ^(C) = and 
goes as follows (see Fig. 9). First, particle q is moved one step North-East {AH{uj) = U). Next, 
also particle p is moved one step North-East {AH{uj) = 3U; see Fig. 9(a)). After that, particle 
p is moved one step South-East {AH{lu) = 2U), and particle q is moved in two steps to the site 
at dual distance 2^2 from particle p in the North-East direction {AH{lu) = 2U: see Fig. 9(b)). 
It is possible to continue following a pattern of this type until particle p is adjacent only to 
the (original) corner particle of type 1 at the end of the bar "just below" p {AH{lo) = 3U] see 
Fig. 9(c)). Call rji the configuration reached after this last step. Particle q can now be moved 
to the site at dual distance 2-\/2 from particle p in the South-East direction (AiJ(cj) = 3U ; see 
Fig. 9(c)). Call this configuration 772. Move particle p first one step South-East {AH(oj) = 3U) 
and then one step South-West {AH{oj) = U). Finally, move particle q to the free site adjacent 
top {AH{uj) = 0). 

Remark 4.4 Note that from 772 to rj' particle p moves in the South direction via the same 
mechanism that was used to move in the East direction from rj to rji. This symmetry in the 
motion of the dimer around a corner of the cluster will be used also in the sequel. 
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(a) (b) (c) 

Figure 9: A dimer is moved to a comer within energy barrier 3U. 



2. The second method is achieved within energy barrier t/ + 4Ai (i.e., max^g^ -^(C) = U + 4:Ai) 
and goes as follows (see Fig. 10). First, move particle q one step in the North-East direction 
(A7?(cj) = [/), and let two extra particles of type 1 enter A and reach the two sites at dual 
distance 1 from q in the West and the South direction {AH{uj) — U + 2Ai). Next, move particle 
p one step in the North-East direction (Ai7(w) = U + 2Ai, see Fig. 10(a)). After that, move 
the particle of type 1 at dual distance 1 in the West direction from p to the site adjacent to 
p in the South- West direction {AH{uj) = U + 2Ai), let one extra particle of type 1 enter A 
{AH{Ld) = 2 + 3Ai), and move this particle to the site at dual distance 1 in the West direction 
from p {AH{uj) = 3Ai, see Fig. 10(b)). 





(d) 



(e) 



(f) 



Figure 10: A dimer is moved to a corner within energy barrier U + 4Ai. The small circles 
represent the extra particles of type 1 with respect to those in the starting configuration. 

Move the particle adjacent to p in the South-East direction one step in the North-East 
direction (AiJ(w) = U -|-3Ai). Movep one step in the South-East direction {AH{bj) = U + 2Ai, 
see Fig. 10(c)). Afterwards, use one of the free particles of type 1 tosaturatep {AH{uj) = +3Ai), 
and remove the other free particles of type 1 from A {AH{u!) = 2Ai, see Fig. 10(d)). The same 
procedure can be repeated until the configuration in Fig. 10(e) is reached {AH{uj) = 3Ai). Next, 
let a particle enter A and reach the site at dual distance 1 in the East direction from p, and 
move one step in the South-East direction the particle adjacent to p in the South-East direction 
(AiJ(w) = U+AAi, see Fig. 10(f)). Next, move p in the Sout-East direction {AH{lj) = U+iAi), 
saturate it with one of the free particles of type 1, and remove the other particles of type 1 from 
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A {AH{lo) = 2Ai). Now particle p can be moved in the South direction in the same way it was 
moved in the East direction within energy barrier U + 3Ai. ■ 



By Lemma 3.3, we know that any path in (□ — > ffl)opt enters the set V*^f*(f*_i)+2 when 
a particle of type 2 is added in d^A to a configuration in g(n, Vj,^^*(£*_i)+i). By Lemma 3.5, 
we know that ^(n, V^,.£*(£*_i)+i) can be determined by looking at all the configurations that 
can be reached starting from the standard configurations in V-i,^e*{e*-i)+i without changing the 
number of particles of type 2 in A and taking into account all the moves that are possible within 
energy barrier A2. Note that different configurations in V^^^.(^*_i)+i necessarily have different 
tile support, and so to move between these classes of configurations it is necessary to perform a 
sequence of heavy-steps. 

Remark 4.5 Note that, from the point of view of the maximal energy barrier that needs to be 
overcome to go from 77 to rj' , for Ai > the first method is more efficient while fOr Ai < 
the second method is more efficient. 

Remark 4.6 Starting from a 2-tiled configuration with a monotone dual support inscribed in a 
rectangle of side lengths li,l2 ("far enough" from the boundary of A), it is possible to reach, via 
one of the two mechanisms described above, all the 2-tiled configurations with a monotone dual 
support and the same circumscribing rectangle. 

4.2 Restriction on heavy-steps 

Lemma 4.7 Let A2 < 31/ + Ai and rj e V^,.£*(£*_i)_|_i. // the first heavy-step starting from rj 
( that does not change the number of particles of type 2 ) does not result in the motion of a comer 
particles of type 2 along the edge where in rj it shares a bond with a comer particle of type 1, 
then it is not possible to reach a new configuration rj' £ V^,.^*(^*_i)+i within energy barrier A2. 




Figure 11: Types of edges for the first possible heavy-step starting from a standard configuration. 

Proof. It is clear that the first heavy-step can only involve one of the external particles of 
type 2. For the proof we refer to Fig. 11. where a prototype configuration 77 e Vi,^^*(f*_i)+i is 
represented. 

We will show that if a heavy-step is performed along one of the edges a, b or 6', then it is not 
possible to reach a new configuration r/' e V^,.^*(^*_i)+i without exceeding energy barrier A2. 
These edges are representatives of the possible types of edges along which the first heavy-step 
is possible without involving the motion of a corner particle of type 2 along the edge where it 
shares a bond with a corner particle of type 1. 

1. Assume that the first heavy-step is along edge a. Let {u,v) denote the dual coordinates of 
the particle pi of type 2 we want to move. Let rji ^ Vi,,£*(£*-i)+i be the configuration that is 
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reached when pi is moved, along edge a, to the site with dual coordinates (u + i,u + i), and 
let ?7o be the configuration visited just before the heavy-step is performed. When tji is reached, 
either the site with dual coordinates (u + 1,?; + 1) is empty or it is occupied by a particle of 
type 1. Without loss of generality, we may assume that r/i does not contain free particles of type 

1, since these particles could be iteratively removed from A while decreasing the energy of the 
configuration. Similarly, we may assume that all particles of type 1 that do not interfere with 
the heavy-step that is performed are still in A, since in fj they had at least one active bond and 
thier removal would increase the energy of the configuration (since Ai < U). In the former case, 
B{r]i) = ^(77) - 5, and so H{r]i) - H{fi) = 5C/ - Ai > A2. In the latter case, B{-qi) = B{f]) - 4, 
and so H{rii) - H{fi) = AU > A2. 

2. Assume that the first heavy-step is along edge 6 (see Fig. 12(a)). Again, let w = {u, v) denote 
the dual coordinates of the particle pi of type 2 we want to move. Denote by qi and (72 the 
particles of type 1 sitting in fj, respectively, at the sites with dual coordinates (u + i, w — i) and 
{u + ■^jV — ■^). Let rji ^ be the configuration obtained by moving particle pi to the 
site X with dual coordinates {u + ~ i), and let rjo be the configuration visited just before 
the heavy-step is performed. When particles pi reaches site x, it has at most two active bonds, 
depending on whether the dual sites yi — (u-|-l, v— 1) and y— 2 = (u+1, v) are empty or occupied 
by a particle of type 1. If both yi and y2 are empty, then Hijji) — H{fj) = hU — 2Ai > A2 (again 
we assume that A does not contain free particles of type 1 and all other particles of type 2 are 
saturated). If only one dual site between yi and y2 is occupied, then, arguing as before, we get 
i/(7/i) - H{fi) = AU > A2 if particle q2 is stiU inside A and H{t]i) - H{f]) = 4t/ - Ai > A2 if 
particle qi has been removed from the A. If both yi and j/2 are occupied, then H{rii)~H(f]) — 3U 
(again we assume that there are no free particles of type 1 in A). Note that, since A2 < 3U -t- Ai, 
no particle of type 1 is allowed to enter A nor is it allowed to break any active bond. Therefore 
the only moves that are possible starting from 771 are those that do not increase the energy. 
Only two moves are possible. Either the particle of type 2 is moved back from site x to site w, 
or the particle q^ of type 1 sitting at the dual site y^ — {u — + is moved to site x. In 
the former case, it is clear that the configuration that is reached is 770 and the move produces a 
non-self-avoiding path. In the latter case, we reach a configuration 772 ^ V*/*(£*_i)+i with the 
same energy as 771 . But now the only move that does not increase the energy is the motion of (73 
back to 773, which again produces a non-self-avoiding path. 



■ ■ ■ ■ : 

■ ■ W ■ 

(a) (b) 
Figure 12: Two choices for the first heavy-step that are too costly. 




3. The case where the first heavy-step is performed along edge b' is similar to the previous case 
(see Fig. 12(b)). ■ 



We see from Lemma 3.6 that ^(n, V*.£*(£*_i)-|_i) is a good candidate for 



bd- 
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4.3 Small values of t 



In Section 5 we will identify the geometry of the protocritical and critical configurations for > 4, 
i.e., for the subregion A2 > 4J7— |Ai. The analysis will show that the set .g(n, V^,^*(f*_i)+i) — 
^(□j V^, consists of all the configurations in Vj,^^*(£*_i)+i. The region A2 < 4J7— |Ai 

is denoted by RT. 

For I* — 2 this set consists of those 2-tiled configurations whose dual tile support is either 
a 2 X 2 square from which a corner has been removed or a 3 x 1 rectangle. For = 3 it 
consists of those 2-tiled configurations whose dual tile support is either a 3 x 2 rectangle plus 
a "protuberance" on one of the four sides or a 3 x 3 square from which two corners have been 
removed. This can be easily verified by noting that it is possible to move a tile protuberance 
within energy barrier AU — 2Ai < A2 and that it is possible to "slide" an external 12-bar of 
length 2 within energy barrier 2U + Ai < A2, as described next (see Fig. 13). 




(a) (b) (c) (d) 



Figure 13: Sliding a 12-bar of length 2. 

We construct a path lj from configuration 77 of Fig. 13(a) to configuration 77' obtained by 
shifting the 12-bar of length 2 in the East direction by one dual unit. Let pi denote the Eastern- 
most particle of type 2 of the 12-bar and p2 the other particle of type 2. Let qi be the particle 
of type 1 adjacent to pi in the North-East direction, (72 the particle of type 1 adjacent to pi 
in the North- West direction, and (73 the particle of type 1 adjacent to P2 in the North- West 
direction. Move qi one step North-East {AH{uj) — U), let an extra particle (74 of type 1 enter A 
{AH{uj) = U + Ai), and let this particle reach the site at dual distance 1 in the West direction 
from qi, and move pi one step North-East (AiJ(aj) = 2U + Ai; see Fig. 13(b)). Then, without 
increasing the energy of the configuration, move p2, qi and q^ subsequently one step South- 
East {AH{uj) = Ai; see Fig. 13(c)). Afterwards, move first (72 and 53 one step North-East 
(Ai?(w) =2U + Ai) and p2 one step North-East (A7f(w) = 2U + Ai; see Fig. 13(d)). Finally, 
move p2 one step South-East, use (72 to saturate p2, and remove (73 from A {AH{uj) = 0). 

It turns out that in each of these cases V — V'*^f*(f*_i)+i and C^^ = V^^?(^,_]^)_|_]^. The proofs 
are essentially analogous to those that will be given in Section 5.1 below. 

5 Proof of Theorems 1.7 and 1.8 1.10 

In Sections 5.1-5.3 we will identify V and C,*^ for the subregions RA, RB and RC, respectively. 
Once the structure of the configurations in V and C^^ are identified, (II3-a) and (II3-b) will follow 
immediately. To prove (H3-c), we will show the existence of a C G C*^^^{fi) such that (!>{(, ffl) < F* 
(i.e., the existence of a "good site" in rj). [CHECK THIS!] Due to the fact that A has a border 
of width 3 where there is no interaction between particles (recall item 3 in Section 1.4), it will be 
immediate that such a good site can always be reached by a particle of type 2 in 9~A without 
touching the cluster. Finally, to see that all a; G (□ — > ffl)opt pass through C*tt, we observe that, 
as long as the free particle of type 2 does not attach itself to the cluster, the energy cannot drop 
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below r* and hence no other particle is allowed to enter A. In particular, this implies that the 
set V*^>f*(f*_i)+2 3 ffl cannot be reached. 

In each of the following sections we look at a specific standard configuration, labelled by 
the position of the lower-left particle, acting as the representative of the set of all standards 
configurations that are obtained by translation and rotation. Each section is split into three 
parts: (1) identification of gdrj}, Vj,^f*(£*_i)+i) and g{{f]}^ Vj,^f*(f*_i)+i) via the motion of 2-tiles 
between dual corners; (2) existence of a good site via the construction of a path to ffl below energy 
level r*; (3) identification of V and C^^. 

5.1 Region RA: proof of Theorem 1.8 

Let Va be the set of 2-tiled configurations with £*{£* — 1) + 1 particles of type 2 whose dual tile 
support is a rectangle of side lengths £*,£* — ! plus a protuberance on one of the four side of the 
rectangle (see Fig. 14). 



Pi- 



ca) 



(b) 



Figure 14: Two examples of configurations in T>a for £* = 7. The external particles of type 2 are 
those lying in the shaded area. 



Definition 5.1 A path u) from 77 € to rj' € V^,^£*(£*_i)+i such that n.2(0 ~ ^^{'n) 

for all ^ € u! is called a modifying path from rj to rj' . 

With this definition, a configuration rj' e V^.f*(f*-i)+i belongs to g{{fj},Vi,^t*{i<--i)+i) (respec- 
tively, g({7?}, V*^£*(£*_i)+i)) with f] a standard configuration in V,t_£*(£*_i)+i if and only if there 
is a modifying path from fj to rf that does not exceed (respectively, stays below) energy level F*. 

Remark 5.2 Note that if there is a path lo: rj ^ rj' that does not exceed (stay below) T* , then 
there is also a path uj' : rj ^ rj' that does the same without ever completing a heavy-step and 
reaching a configuration that is equivalent to a configuration that has already been visited. 

5.1.1 Identification of g{{fj},V^j^e-^i)+i) and g({?7}, 
Lemma 5.3 .g({^}, V*/*(£*-i)+i) = g({?7}, K/*(i^*-i)+i) ='Da ■ 

Proof. Let p be a configuration consisting of 2-tiled dual rectangle with horizontal side length 
£* and vertical side length £* — 1 whose top-rightmost 2-tile is centered at site (op, bp). Let fj be 
the 2-tiled standard configuration obtained from p by adding a 2-tiled protuberance centered 
at (dual) site x — [op + 1,6^ + 1) (see Fig. 11). Observe that H{fj) = F* — A2. Let y be 
one of the two sites of the tile centered at x that is occupied by a particle of type 1 with only 
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one active bond. It is easy to check that Va C 5({77}7 Indeed, it is enough to 

consider the configuration (p, y) obtained from fj by first detaching and removing the two corner 
particles of type 1 adjacent to x and then moving the particle of type 2 from x to y, which gives 
^{{p,y),ri) < r* for all rj G Va- Note that this is true irrespective of the distance of ?/ to the 
boundary of A. 

To conclude the proof we will show that all modifying paths starting from fj either lead to a 
configuration in Va or exceed energy level T* . Note that we are interested only in those modifying 
paths consisting of at least one heavy-step, otherwise the only configuration in Vj,^*(f*_i)_|.i that 
can be reached from ^ is 77 itself. In other words, since all configurations in Va consist of a 2-tiled 
dual rectangle plus a 2-tiled protuberance, we have to show that, without exceeding energy level 
r*, i.e., within energy barrier A2, the only configurations in Vj,_£*(£*_i)+i that can be reached 
via a path consisting of configurations with — 1) + 1 particles of type 2 are obtained by 

"moving" the 2~tiled protuberance. 

Let rj G 3({?7}, V,t_f*(£*_i)+i); and let w be a modifying path from fj to rj such that H{^) < 
r* H{7j} + A2 for aU ^ e uj. Note that, since A2 < 3U, we have B-{^) < 2 for all ^ e uj. 
Furthermore, in VaMVa n RT) we have Ai > and A2 > ?7 + 3Ai. This imphes that, for 

(i) if B-{^) = 1, then < 2; 

(ii) if B- {1) = 2, then < 1. 

Fig. 11 shows the different classes of edges along which the first heavy-step of a modifying path 
is possible. We will group modifying paths according to the edge along which the first heavy-step 
is made. Note that A2 < 3U + Ai in region RA and so Lemma 4.7 applies. Therefore, as a first 
possible heavy-step, we need to consider only those consisting of the motion of a corner particle 
of type 2 along the edge where it shares a bond with a corner particle of type 1. In order to 
identify particles and edges, we refer to Fig. 11. 



■ : • : ■ : • 




(a) (b) (c) 



Figure 15: Region RA, first heavy-step along edge c. 

Step 1: 

Claim 5.4 All modifying paths starting with a heavy-step involving a particle of type 2 other 
than that in the protuberance exceed energy level T* . 

Proof. 1. Asumme that he first heavy-step is along edge c. Let x be the site where the particle 
pi that we want to move sits in configuration fj. Let ?7i_ be the configuration visited just before 
the first heavy-step is performed, and let 771 be the configuration that is reached when pi is 
moved one step North-East to site y. There are four possible cases. 

• B(pi,rii) — 0. Let us assume that all free particles of type 1 are removed from A. Then 
H{rii) — H{fj) = — Ai > A2 and hence the path exceeds energy level F"^. 

• S(pi,77i) = l: i/(r/i)-i/(?7)=3C/> A2. 
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• -B(pi,r/i) — 2. Let us assume that the two sites occupied by a particle of type 1 are y\ 
(South-East of y) and j/2 (North-East of y). Since Ai < U ^ the least expensive way in terms 
of energy cost to have B(j)\,r\\) — 2 is achieved by bringing one extra particle of type 1 
in A (see Fig. 15(a)): H{r]i) — H{'q) = 2U + Ai. Therefore, starting from r^i, only moves 
that do not increase the energy are possible. Since 771 is not equivalent to any configuration 
in Vi,^£*(£*_i)+i, at least one extra heavy-step is necessary from 771. By Remark 5.2, the 
first heavy-step from rji cannot be completed by moving pi back to x. Before the next 
heavy-step is performed, the only moves from rji that do not increase the energy further 
are motions of particles of type 1 with one active bond to or from a site adjacent to pi. In 
particular, it is not possible to bring inside A any other particle or type 1. Any possible 
sequence of such moves cannot change the energy of the configuration. When the next 
heavy-step is completed, at least one extra bond is added and energy level T* is exceeded. 

• B{pi,rji) = 3 (see Fig. 15(b)). Similarly to the previous case, the least expensive way to 
have _B(pi,7/i) = 3 is achieved by bringing two extra particles of type 1 inside A. From 
771 at least one other heavy-step is necessary, and moves that increase the energy are not 
allowed. As in the previous case, only motions of particles of type 1 that do not decrease 
the number of bonds are allowed (see, for instance. Fig. 15(c)), and the completion of the 
next heavy-step exceeds energy level F"^. 




(a) (b) (c) 

Figure 16: Region RA, first heavy-step along edge d. 



2. Assume that the first heavy-step is along edge d. Let x he the site where the particle pi that 
we want to move sits in configuration 77. Let ?7i_ be the configuration visited just before the 
first heavy-step is performed, and let rji the configuration that is reached from 77 after moving 
pi one step North- West along the edge d to site y. Also, let p2 denote the particle of type 2 in 
the tile protuberance, and p^ the particle of type 2 at dual distance one from pi in the South 
direction in configuration fj. As in the previous case, after the heavy-step is performed we must 
have B{pi, 771) > 2. 

• B{pi,r]i) = 3. The best choice is when B^ [rji) ~ 1 and nl{rii) — 2 (see Fig. 16(a)). From 
771 only moves that do not increase the energy are possible. There is only possible one non- 
backtracking move: the particle of type 1 South- West of x is moved one step North-East. 
After this move, it is only possible to move the particle of type 1 South- West of pi one 
step South- West. The configuration that is reached does not belong to V*.^*(^*-i)+i and 
no other move is allowed (this is the analogue of i3(pi, 771) = 3 in the previous case). 

• B{pi, rji) = 2. As in the previous case, the choice that minimizes the energy of 771 is such 
that B^ [rji) — 2 and 7^^(771) = 1, and the particles of type 1 adjacent to pi sit at sites yi 
and ?/2, respectively, South-West and North-West of y (see Fig. 16(b)). It is easy to see 
that this is the choice that allows for "more freedom" of the path, in the sense that it is 
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the only choice from which it is possible to complete a further heavy-step. Therefore, only 
moves that do not increase the energy are allowed starting from rji, and again note that 
no other particle of type 1 is allowed to enter A. From 771 at least one other heavy-step is 
necessary. This means that from rji it is only possible to move pi one step North-East to 
site z (this case will be examined afterwards) or to start a sequence of moves of particles 
of type 1 to or from a site adjacent to pi. Any possible sequence of such moves cannot 
decrease the energy of the configuration and only a heavy-step can be completed by moving 
Pi one step North-East to site z to reach configuration 772. Clearly, the configurations the 
path can reach (strictly) below F* from 772 are the same as those that can be reached from 
773 by saturating pi with the two free particle of type 1 (see Fig. 16(c)). Arguing as in the 
previous case, we see that the only heavy-step possible without exceeding energy level F"^ 
from 773 is the one that is completed by moving pi back to site y (configuration 774). The 
transition from 773 to 774 can be treated in the same way as that from 77 to 771, and so we 
conclude that from 774 it is only possible to reach a configuration equivalent to fj. 



Step 2: 

Claim 5.5 Let the first step of a modifying path starting from a configuration in Va involve the 
particle p of type 2 in the protuberance. If p is not re-attached to the main cluster before the next 
heavy-step is completed, then the path exceeds energy level T* before reaching a configuration in 

Proof. The proof of this claim will be deferred to the proof of Lemma 5.9, Step 1, in Section 5.2. 
There it will be shown that the same claim holds also in region RB, where the values that A2 
can take are larger than those that A2 can take in RA. ■ 

Depending on its starting position, the protuberance can be re-attached in two possible ways: 
either the particle of type 2 shares two particles of type 1 with the other particles of type 2 in the 
cluster (see Fig. 17(a)), or it shares only one particle of type 1 (see Fig. 17(b)). In both cases we 
consider the evolution of the path to a 2-tiled configuration that is equivalent to the one reached 
the moment the particle of type 2 joins the main cluster. In the first case, the configuration is 
again in the class Va^ and hence the same kind of argument can be repeated. In the second case, 
the following statement allow us to conclude the proof of Lemma 5.3. 



(a) (b) 

Figure 17: Modifying paths may create hanging protuberances. 

In the first case, the configuration is again in the class Va and, hence, the same kind of 
analysis can be repeated. In the second case, the following statement allow us to conclude the 
proof of the lemma. 

Step 3: 
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Claim 5.6 From a configuration consisting of a 2-tiled rectangle plus a hanging protuberance it 
is not possible to reach a configuration in without exceeding energy level T*. 

Proof. Let rjo be the configuration of Fig. 17(b), r]i the configuration reached when the first 
heavy-step from r]o is completed, and r]i- the configuration visited by the path just before r]i. 
Let pi be the particle of type 2 in the hanging protuberance and p2 the particle of type 2 at dual 
distance in the South-East direction from pi. Let gi denote the particle of type 1 shared by 
pi and P2- We will show that if the first heavy-step from r]o is not completed by moving pi one 
step North-East of one step South- West (the two cases arc analogous), then the path exceeds 
energy level T*. From Lemma 4.7 and Step 1 it follows that, in order to prove the claim, we 
need to consider only the heavy-steps completed by moving pi North- West or South-East and 
P2 North- West. 

• Assume that pi is moved one step South-East. Observe that B{pi,r]i) < 2 and B{p2,r]i) < 
3. Clearly, since Ai < U, the choice that is most favorable from the point of view of 
energy is when pi has 2 active bonds and p2 has 3 active bonds. It is clear that, since rji 
is reached from rji- via the motion of a particle of type 2, the two configurations have the 
same particles of type 1 placed at the same sites. Since B{pi,rii) = 2, there is no advantage 
in having more than two particles of type 1 adjacent to pi in rji-, since one bond will be 
lost anyway with the motion of p2. It follows that H{r]i) — H{f]) = 3U + Ai > A2 (see 
Fig. 18(a)). Arguing in the same way, we see that in the case where the first heavy-step from 
rjo is completed by moving p2 one step North- West, we have H(ri\) — H{ff) = 3?7-|-2Ai > A2 
(see Fig. 18(b)). (Note that 3f7 -|- Ai > A2 also in region RB, and so these moves will be 
forbidden there as well.) 

• Assume that the first heavy-step is in the North- West direction. Then the only possibility 
without exceeding energy level F* is when B{px,r]\) = 2. This is achieved, for instance, 
by moving the two particles of type 1 on the West side of p\ one step North- West before 
completing the heavy-step. Since 771 has two broken bonds, it is not possible to break any 
extra bond, and hence no extra heavy-step is possible as long as the particle of type 1 does 
not reach a site adjacent to p\. Let r/2 be such a configuration (see Fig. 18(d)). Since r/2 
has one broken bond and one extra particles of type 1, the next heavy-step (completed by 
moving either pi or P2) cannot break more than one bond, but clearly this is impossible. 

• Assume that the first heavy-step is completed by moving pi North-East (see Fig. 18(c)). 
Then rji is a configuration that can be reached with one heavy-step from a configuration in 
"Da; and hence the claim in Step 2 holds. This can be done by moving the two particles of 
type 1 North of px one step North-East (breaking two active bonds) , moving pi one step 
North-East {AH = 0) and removing the free particle of type 1 froma A {AH = — Ai). 



5.1.2 Existence of a good site 

Lemma 5.7 For all r) e Va, there exists an x & F{fj) such that ??((f},x),ffl) < F*. 

Proof. Configurations in Va consist of a dual 2-tiled square and a protuberance on one of the 
four sides. If the protuberance is on the longest side, then the rectangle circumscribing the cluster 



28 



(a) (b) (c) (d) 

Figure 18: Region RA, first heavy-step along edge e. 



is a square of side length £*. Conversely, if the protuberance is on one of the shortest sides, then 
the rectangle circumscribing the cluster has side lengths i* — + 1. Without loss of generality, 
we assume that the longest side of the rectangle is horizontal and that the protuberance is on 
the North side of the rectangle when it is attached to the longest side and on the East side when 
it is attached to the shortest side. 

• Assume that the protuberance is on the longest side. Let x be the central site of a tile 
adjacent to the protuberance to the North side of the rectangle. After a particle of type 2 
has entered A and has been moved to x, a configuration with energy T* ~3U is reached. The 
particle at site x can be saturated within the energy barrier Ai, to reach a configuration 
with energy T* — A2 — s < H{ri) for all rj € Va, consisting of a rectangle of side lengths 
£*,£* — ! plus a 12-bar of length 2 on its (longest) north-side. Then, within energy barrier 
A2, the 12-bar can be completed to obtain a 2-tiled dual square of side length £*. The 
claim in the lemma follows from Remark 2.1. Such a configuration is supercritical (and 
has energy smaller than the standard configuration with £*{£* ~ 1) particles of type 2). 
Therefore it can grow until ffl is reached while staying below energy level T* . 

• Assume that the protuberance is on the shortest side. Let x be the central site of a tile 
adjacent to the protuberance to the East side of the rectangle. After a particle of type 2 
has entered A and has been moved to x, a configuration with energy F* — 3U is reached. 
The particle at site x can be saturated within energy barrier Ai, to reach a configuration 
with energy T* — A2 — e < H{rj) for all rj G Va, consisting of a rectangle of side lengths 
£*,£* — 1 plus a 12-bar of length 2 on its shortest (see Fig. 19(a)). 




(a) (b) (c) 

Figure 19: Region RA, first heavy-step along edge e. 



Let y be the central site of a tile adjacent to the North side of the rectangle of side lengths 
£*,£* — !. One of the dimers of the 1 2-bar of length 2 can be moved within energy barrier 3C/ 
(note that 3C/ — A2— e < 0) to the north Side of the rectangle (as described in Section 4.1), 
with the particle of type 2 at site y, to reach a configuration with energy T* — A2 — e + U. 
The particle of type 2 can be saturated using one of the corner particle of type 1 adjacent 
to the protuberance that is left on the East side. This is done within energy barrier U, and 
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the resulting configuration again has energy F* — A2 — £ + [/ (see Fig. 19(b)). Finally, move 
the pending dimer that is left on the East side to a corner adjacent to the tile centered at 
site y (see Fig. 19(c)). As described in Section 4.1, this can be done within energy barrier 
2U (since the particle of type 2 of the dimer has only one neighboring corner particle of type 
1). The configuration obtained in this way consists of a rectangle of side lengths — 1 
plus a 12-bar of length 2 on its longest North-side. It is now possible to continue as in the 
previous case. Note that, to move the 12-bar of length 2 from the East side to the North 
side, we implicitly assumed that these sides of the rectangle are far from the boundary of 
A. 

Now suppose that the protuberance was originally on a side of the rectangle close to A. 
Then we can proceed in the following way. After the 12 bar of length 2 has been completed, 
it is possible to complete the 12 bar below energy level F* and obtain a 2-tiled rectangle 
of side lengths £* + — I. From this configuration it is possible to iteratively remove 
— 3 corner 2-tiles from the shortliest side of the rectangle that is far from A, in order 
to obtain a 2-tiled configuration where the 12-bar of length 2 is far from d~A as well. 
Clearly, this can be achieved below energy level F*. It is now possible to proceed as in the 
previous case, observing that either the North side or the South side of the rectangle is far 
from d~A. 



5.1.3 Identification of V and C^^ 
Lemma 5.8 V =Va and C^^ = V^^"^. 

Proof. All configtirations in have energy F*. We will show that all configurations in T>^^ 
are essential saddles, and therefore belong to G{0, ffl) by Lemma 1.3. Pick 77 — {fj, x) in V^"^. 

Let U{r]) be the set of optimal paths entering V^^^*(^*_i)_|_2 via configuration rj. Pick w €U{r]) 
such that H{^) < T* for aU ^ e P^v) and S^{ri) = {{fj, yi), ...,{ri, ym), • • • , H}, with yt ^ d'A 
for all j G 1. • . • ! m and in F{fj) such that H{^) < F* for all ^ G Su;{{fj,ym))- Note that such 

a path Lo exists because: 

(i) f] e g(n, V^_£*(£,_i)+i) by Lemma 5.3; 

(ii) there exists an y^ € F{ri) such that <P((fj, y,n)) < F* by Lennna 5.7; 

(iii) in configmration 77 there exists a lattice path of empty sites yi,...,ym with yi x and 
yi ^ d~A for alH € 1, . . . , m, since for all x G d~A there exists an y ~ a; such that y ^ d~A 
(by the shape of A), and [fj] is "sufficiently far" from d~A (i.e., there is a border of width 

3 in A where particles do no interact). 

Next, note that all saddles in Suj{rj) are configurations of the type 77' = {fi,u), u ^ d~ A. Recall 
Definition 1.2(e). We have to show that S'(w') ^ S{ijo)\{r}} for all w' € (??, ??')opt- Consider 
the partition (□ ffl)opt = (U(ri),W{ri)) with W{ii) = (□ ^ ffl)optW(?7). If uj' G U{r]), then 
S{uj') ^ S{uj)\{r]} because 77 G uj' . If, on the other hand, w' G W^ir]), then, by Lemma 3.3(1), ui' 
enters the set Vi/*(f*_i)+2 via some configuration ( = {C,z) with H{() = F* and z G d~A that 
does not belong to S(uj) (by the construction of ui). Since the choice of 77 G V^'^ was arbitrary, 
we conclude that all configurations in are essential saddles, and hence that C G{0, ffl). 

To prove that V^'^ — C^^, i.e., is the entrance of the essential gate, we will show that, 

for any path cj G (□ ^ ffl)optj any configuration ( G S{0, ffl) that is visited by lo before some 
configuration in V^"^ is an unessential saddle and therefore does not belong to G{0, ffl). 
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• We show that ah saddles in V*_<f*(f*_i)+i are unessentiaL To that end, we pick S (□ — 
ffl)opt and we let C e w be such that 7?(C) =r* and 713(0 <i*{e*-l) + l- By Lemma 3.3(1), 
all optimal paths enter V^.£*(^*_ij+2 via a configuration in g{0, V*/*(f*-i)+2)'^'^^- Hence, af- 
ter visiting (, uj must visit g{0, V*,£*(£*_i)+2) before entering V^_£*(£*_i)+2- By Lemma 5.3, 
combined with Lemma 3.5, we have ^(n, Vi,f*(f*_i)+2) = g(,OT^*,e*{e*-i)+2)- The claim 
now follows via Lemma 3.6(1). 

• For all w e (□ — > ffl)opt such that ^ £ S{uj) and ^ G V'^,>£*(j?*-i)+2\2-'l4'^^i there exists an 
T] G such that 77 G Plj{0- Indeed, by Lemmas 3.3(1) and 5.3, all optimal paths enter 
V*,£*(£*-i)+i via a configuration in T)^^. 

• To conclude, we need to show that V — 'Da- The inclusion V ^ T^a is immediate, since 
all rj — (r), a;) G T)^^ = C^^ are obtained, in a single step, by adding a particle of type 
2 at the boundary of A to a configuration 77 G (?(□, Vi,_£*(£*_i)_|_i) = Va- To see that 
V C Va, suppose that there is a configuration 77 G VXD a- Let fj' be the configuration in 

= I?^*^^ obtained in a single step from fj. Observe that C r^V^ 1)^2 • Since, by 
Lemma 3.3(1), all optimal paths enter V^.£*(f*_i)+2 by adding a particle of type 2 at 
to a configuration in g(n, Vj,^£*(£*_i)-|-i) — T^a, it follows that f] G V*^>f*(f«_i)+2. Indeed, 
by assumption, this configuration cannot be in (7(n, Vj,^£*(£*_i)+i) = T^a- In particular, it 
follows that f] G V*^>f*(f*_i)-)_2\2-'A^^- But this means that there is a path w G (□ — > ffl)opt 
that reaches ?) G V*_>f*(f*_i)+2\2-'j4'^^ and does not contain any configuration in t/(n,ffl). 
In particular, Cj does not contain any configuration in T)^^. But this is a contradiction, 
since by Lemmas 3.3(1) and 5.3, all optimal paths enter via a configuration 

in I?bd2 



5.2 Region RB: proof of Theorem 1.9 

Let I?B be the set of 2-tiled configurations with £*{£* — 1) + 1 particles of type 2 whose dual 
tile support is a monotone polyomino and whose circumscribing rectangle has side lengths either 
£*,£* or r + l,r - 1 (see Fig. 20). Note that Vb 2 Va- 




(a) (b) 

Figure 20: Two examples of configurations in Vb for £* = 7. 



5.2.1 Identification of g{{v},V^j*{e*-i)+i) and g{{v},V*,e*{e*-i)+i) 

There are configurations in Vb that cannot be reached within energy level T*. These config- 
urations have support near the boundary of A and for A — > form a negligible fraction of 
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Vb. 

Lemma 5.9 V*/*(£*-i)+i) = V*/*(f*-i)+i) Q Vb- 

Proof. Note that in P^, for ^ G Vj,_^*(£*_i)+i the fohowmg conditions are satisfied 

(i) B-{<j) = l: n+{a)<3; 

(ii) B-{(t) = 2: n+(cr) < 1. 

(iii) B-(ct) = 3: n+{a) < 0. 

Observe that A2 < 3U + Ai throughout region RB, and therefore Lemma 4.7 apphes. This 
means that any heavy-step, completed without exceeding energy levl T* and starting from a 
configuration in V^,^£*(£*_i)+i, necessarily involves a corner particle of type 2 that is moved along 
the edge where it shares a bond with a corner particle of type 1. 

From Remark 4.6 it follows that all 2-tiled configurations with a monotone dual support, 
a circumscribed rectangle of side lengths {£*,£*) or {£* ~ 1,£* + 1), and with a fixed lower left 
corner far enough from d~A, belong to (7(f?, V*,f*(f*_i)+i)- It remains to show that all other 
configurations in V*_f*(£*_i)+i cannot be reached without exceeding energy level F*. As in the 
case of the analogous lemma for region RA, the proofs comes in various steps. 

Step 1: Let the first step of a modifying path starting from a configuration in Va involve the 
particle p of type 2 in the protuberance. If p is not re-attached to the main cluster before the next 
heavy-step is completed, then the path exceeds energy level F* before reaching a configuration in 
Vi,^f*(f*_i)_(_i. Therefore we may consider a path whose first heavy-step involves the protuberance 
and whose second heavy-step does not. 




(a) jyi (b) ri2 (c) r/a (d) r;4 (e) r/5 



Figure 21: The presence of a detached 2-tile precludes the motion of other 2-tiles. 

Let be a standard configuration whose protuberance on the North side belongs to the 
Western-most bar, and let 771 (see Fig. 21(a)) be the configuration reached by completing a 
heavy-step that moves the particle pi in the protuberance one step North- West. Since the second 
heavy-step of the path does not involve a motion of pi , we may assume that from configuration 
772 (see Fig. 21(b)) the path saturates particle pi with two extra particles of type 1. It will 
become clear later on that this choice for 770 and rji is the most interesting, since after the next 
heavy-step is completed the particle of type 2 that is moved can share two particles of type 1 
with pi. 

Let p2 be the North- West particle of type 2 in the 2-tiled rectangle of 772, Ps the particle of 
type 2 below p2 ■ and 754 the particle of type 2 East of p2 . It is easy to check that any heavy-step 
completed by moving any other corner particle of type 2 other than p2 leads to an energy above 
F*. Since 772 has two extra particle of type 1, the next heavy-step must be completed by breaking 
at most one bond. (Note that in 772 all external particles of type 1 are saturated and hence at 
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least one bond must be broken.) This means that the next particle of type 2 must be moved 
from one good dual corner to another, and the corner must be created by breaking exactly one 
bond. It follows that the only heavy-step from 772 that can be completed is the one obtained by 
moving p2 one step North- West with a particle of type 1 sitting South- West of it (configuration 
773; see Fig. 21(c)). From 773 it is not allowed to break other bonds. As a consequence, the only 
heavy-steps that are possible are completed by moving p4 South-East or p3 North-East. These 
two cases are analogous. We describe the second one. 

When p2 is moved North-East, the total number of bonds cannot be decreased. This can 
be achieved by moving one step North- West the particle of type 1 adjacent to p3 in r/3, and 
with the help of one extra particle of type 1 reaching the South- West of the "destination" of p3 
(configuration 774; see Fig. 21(d)). From 774 it is not allowed to increase the energy. Since in 774 
all particles of type 2 have at least three active bonds, the motion of the next particle of type 2 
must be from a good dual to another. But from 7/4 only the motion of corner particles of type 1 
is allowed (configuration 775; see Fig. 21(e)), and it is not possible to create a good dual corner 
without bringing a further extra particles of type 1 inside A. Hence the path cannot be extended 
with another heavy-step without exceeding energy level F*. 

Step 2: In this step we will consider the evolution of those paths that visit a 2-tiled configuration 
having a hanging protuberance. For definiteness, we take 77 to be a configuration consisting of 
a 2-tiled rectangle plus a hanging protuberance at the North- West corner (see Fig. 17(b)). Let 
Pi be the particle of type 2 in the hanging protuberance, and let p2 be the particle of type 2 
South-East of pi. 

If the first heavy-step from r] is not completed by moving pi in the North- West, South- West 
or North-East direction, then the path exceeds energy level F* before reaching a configuration 
in V*^£«(£*_i)+i. From Step 3 in the proof of Lemma 5.8, we already know that it is not possible 
to move pi South-East nor p2 North- West. We will investigate what happens when some other 
corner particle of type 2 is moved to complete the first heavy-step from rj. 

Suppose that the first heavy-step is completed by moving particle p3 of type 2 in the North- 
East corner tile of the rectangle (configuration 771). It is easy to see that this move is only possible 
below energy level F* when, in t^i, p^ has three adjacent particles of type 1. This can be achieved 
by bringing inside A two extra particles of type 1 (see Fig. 22). From 771 it is not possible to 
further increase the energy without exceeding energly level F*, and it is immediate that no other 
heavy-step is allowed. 




Figure 22: The presence of a hanging protuberance precludes the motion of other 2-tiles. 

Step 1 and 2 imply the following. A modifying path that does not exceed energy level F"^ 
and starts with any heavy-step involving the particle of the protuberance must "go back" to a 
configuration equivalent to a configuration in Da before some other particle of type 2 can be 
moved. 

Step 3: Let the first heavy-step from a configuration 7^0 G 2?_b be completed by moving a corner 
particle p of type 1 belonging to a 12-bar of dual length m > i. If p is not re-attached to the 
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main cluster before the next heavy-step is completed, then the path exceeds energy level T* 
before reaching a configuration in Vj,.£*(^*_i)+i. 




(a) rio (b) 'iji (c) 772 (d) 773 



Figure 23: First heavy-step from a bar of length > 3. 

Let r/o be as in Fig. 23(a), and let pi, p2 and be the first three particles of type 2, starting 
from the West, of the Northern-most 12-bar. Let the first heavy-step from r/o be completed 
by moving pi North- West. Since we are assuming that the second heavy-step is not completed 
by moving pi, we may consider the path from rji that is obtained by saturating pi with two 
extra particles of type 1 and by moving the particle North- West of p2 one step South- West (see 
Fig. 23(b)). Note that iji can be reached within energy barrier U + 3Ai. 

Since from rji it is not allowed to break another extra bond, the only possible heavy-step is 
the one obtained by moving p2 North- West after an extra particle of type 1 has entered A and 
has reached a site adjacent to the destination of p2 (configuration 772; see Fig. 23(c)). From 772 it 
is not possible to further increase the energy without exceeding level F*, and a further heavy-step 
is therefore not possible. 

Remark: Note that the key observations here are the following. After pi has moved, two 
extra particles of type 1 rae required to saturate it. The presence of two extra particles of type 
1 forces the path to evolve without breaking extra bonds. Since all the external particles of 
type 2 other than pi have at least three active bonds, the motion must necessarily be towards a 
good dual corner (a site with three neighbors occupied by a particles of type 1). Note that when 
a particle is moved, it changes its parity and only particles with the same parity can interact. 
When the second particle p2 of type 2 is moved, it takes the parity of the particles in the tile pi 
belongs to. But after p2 is moved, it can share at most two particles of type 1 with pi, and hence 
another particle of type 1 is neeede. This leads to a configuration with 3 extra particles of type 
1 and one broken bond, and hence it is not allowed to make moves that increase the energy. In 
particular, it is not allowed to bring inside A other particles of type 1, and only particles of type 
1 with one active bond can be moved. Note that, again, all particles of type 2 have at least three 
active bonds. Hence, a further heavy-step would only be possible if a good dual corner can be 
created close to a particle of type 2 without decreasing the number of bonds. These observations 
are key in order to explore what configurations can be reached by a modifying paths without 
exceeding F*. 

Step 4: Let the first heavy-step from a configuration 770 G "Db be completed by moving a corner 
particle pi of type 2 belonging to a 12-bar of dual length 7n = 2. If the path does not exceed 
energy level F*, then one of the following must happen: 

• Particle p is re-attached to the cluster before any other particle of type 2 is moved. 

• If pi denotes the other particle of type 2 in the bar, then the path reaches a configuration 
r]d where pi and p2 are saturated (with the help of three extra particles of type 2), belong 
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to the same cluster and are at dual distance -\/2 from its location in ryo (see Fig. 24(c)). 
From r]d-, the first heavy-step (without backtracking) must be completed by re- attaching pi 
to the cluster. The configuration that is reached in this way is again a configuration that 
can be reached with a single heavy-step completed by moving a corner particle of type 2 
belonging to a 12-bar of dual length m = 2 (see Fig. 24(d)). 




(a) rio (b) rii (c) rja (d) »?2 



Figure 24: First heavy-step from a bar of length 2. 

Let pi be the Western- most particle of type 2 in the bar of length 2 in configuration 770, 
and p2 the particle of type 1 East of pi (see Fig. 24(a)). Let the first heavy-step from rjo be 
completed by moving pi North- West. Suppose that pi is not re- attached to the cluster before 
the next heavy-step is completed. As in Step 3, we can consider the path from 771 obtained by 
saturating pi with two extra particles of type 1 and by moving the particle North- West of p2 one 
step South- West (see Fig. 24(b)). Note that 771 can be reached within energy barrier U + 3Ai. 

Since from t]i it is not allowed to break another bond, the only possible heavy-step is the 
one obtained by moving p2 North- West after an extra particle of type 1 has entered A and has 
reached a site adjacent to the destination of p2- After that, since it is not possible to further 
increase the energy, before the next heavy-step is completed it is necessary to move a particle 
of type 1 to the empty site adjacent to p2, reaching configuration rjd- The first heavy-step from 
r]d, since it contains three extra particles of type 1, must be completed breaking at most one 
bond. This means that a particle of type 1 must be moved from a good dual corner to another. 
This is only possible, without backtracking, by moving pi one step South- West to site x after 
the particle of type 1 sitting at x in T]d has been moved one step North- West. Particle pi can be 
saturated with a free particle of type 1, to reach configuration 772 of Fig. 24(d). Note that 772 is 
the "mirror image" of rji, and hence the same arguments can be repeated. 

For further reference, let us consider also a configuration with a vertical 12-bar of length 2. 
Such a configuration rj^ does not belong to Vb , but could be reached by a modifying path below 
energy level F* (see Fig. 25(a)). 




(a) 77„ (b) (c) r]a: (d) rjy 



Figure 25: Protuberances attached to 12-bars can be treated as 12-bars of length 2. 

The analysis is completely analogous to the case of a horizontal 12-bar of length 2, and the 
claim is that if the first heavy-step is completed by moving the particle pi of type 2 on top of the 
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vertical 12-bar and the path does not exceed the energy level F*, then either pi is re-attached to 
the cluster before the next heavy-step is completed, or after pi has been saturated (configuration 
?7m; see Fig. 25(b)) the next heavy-step must be completed by moving the other particle p2 of 
type 2 originally in the 12-bar to reach a 2-tiled configuration rj^ with three extra particles of 
type 1 where the 12-bar is "floating" on a side of the cluster (see Fig. 25(c)). From t]x, the 
only heavy-step that can be completed is the one obtained by re-attaching p2 to the cluster. 
Then p2 can be saturated to obtain configuration rjy (see Fig. 25(d)). Note that rjy is the mirror 
image of 77^, and the same argument can be repeated. Note that when pi is re-attached, either a 
configuration with again a vertical 12-bar of length 2 is reached, or the path visits a configuration 
where pi, after being saturated, belongs to a hanging protuberance. We will see below that in 
this case the next heavy-step must be completed by moving pi again. 

Thus, after the first heavy-step from a configuration r/o in 2?b is completed by moving a 
particle pi of type 2, this particle must be re-attached to the cluster. This particle can be 
saturated to reach a configuration rji in which pi belongs to a (possibly hanging) protuberance. 
Note that 771 contains one extra particle of type 1, and therefore H{r]i) = H{r]Q) + Ai. Note that 
in the case pi belonged to a bar of length 2, and configuration rji contains two protuberances. 

We claim that, from rji, the first heavy-step must be completed by moving a particle of type 
1 in one of the protuberances. This can be seen as follows. 



■ • U U U 




• - ■- ■ 



(a) rji (b) ri2 (c) r/i. 

Figure 26: Heavy-step from a configuration with protuberances non in TJa 

1. Let rji be a configuration like the one in Fig. 26(a), and let the first heavy-step from 771 be 
completed by moving the particle p2 of type 2 sitting at the Northern-most site of the West bar 
of 771 one step North- West. Let 772 be the configuration that is reached when this heavy-step is 
completed (see Fig. 26(b)). Since r/p has already one extra particle of type 1, it is not possible 
to have B{p2; 772) = 1 (i.e., three broken bonds). Similarly, if we consider the case B{p2, 772) = 2, 
then we must require one extra particle of type 1, but this is incompatible with a configuration 
with two broken bonds. Therefore the only possibility is B(j)2, 7/2) — 1, where the three bonds of 
pi in 772 are obtained with the help of two extra particles of type 1 . From 772 it is not possible to 
increase the energy further. Since in 772 all configurations have at least three active bonds, the 
next heavy-step must be completed by moving a particle of type 2 from a good dual corner to 
another without increasing the number of broken bonds. But this is impossible from configuration 

2. Clearly the same conclusion can be reached also if pi belongs to a hanging protuberance (see 
Step 2). Note that, by completing a heavy-step from 7^1 by moving pi, the configuration that is 
reached is either one heavy-step away from being attached in a good dual corner to form a bar of 
length 777 > 2 (and we know that it must be re- attached by Step 3 and 4), or it is far from good 
dual corners. In this case, if we assume that the next heavy-step is not completed by moving 
again pi , then we can consider the path from the configuration 7/2 obtained by saturating pi (see 
Fig. 26(c)). This configuration has three extra particles of type 1, and hence no heavy-step can 
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be completed without exceeding energy level F* by moving a particle of type 2 different from pi . 
It is straightforward to see in the next heavy-step pi must indeed be re-attached to the cluster. 

Step 5: It follows from Steps 1-4 that the set of single cluster configurations (and, consequently, 
the set of configurations in V^.^*(^*_i)+i) that can be visited by a modifying path that does not 
exceed energy level T* coincides with the set of configurations that can be reached by a modifying 
path whose configurations do not have more than one particle of type 2 not belonging to the main 
cluster such that if a configuration has a particle of type 2 that is not connected to the cluster, 
then this particle is re-attached to the cluster in the next heavy-step, and afterwards is saturated. 
This means that the configurations in Vj,^£*(£*_i)+i that are reached without exceeding energy 
level r* are obtained by iteratively moving a corner 2-tile around the cluster (with the help of 
one or two extra particles of type 1 used to saturate the particle of type 2). We already saw in 
Steps 1-2 that a modifying path starting from a configuration in Va with a heavy-step involving 
the particles of type 2 in the protuberance cannot leave the set Va- In the other cases, if the 
2-tile that is moved reaches a corner, then a new configuration in T>b is reached. Otherwise, a 
configuration with one or two protuberances not belonging to Va is reached. But in this case, 
the path can only proceed by moving one of the protuberances, which eventually reach a corner 
and again producing a configuration in 

Step 6: From what has been seen so far it follows that the set of single cluster 2-tiled con- 
figurations that can be visited by a modifying path starting from a standard configuration fj 
without exceeding energy level F* consists of those configurations that can be reached by either 
moving the protuberance of a configuration in Va, or by iteratively moving corner 2-tile to some 
other corner possibly created by adding an extra particle of type 1. This observation implies 
that g{D,Vi,^i*(i*-i)+i) C T>b- Furthermore, if ?7 e Vb and has support far enough from d~A, 
then rj G g(n, £*(£*_i)+i). It is straightforward to see that lattice distance 2 from the annulus 
where no interaction is present is already far enough. In order to prove that g{0, V* = 
g{0,%,e-{e*-i)+i), we will show that the set g{0,%^e*{e*-i)+i)\g{0,%,e*{e*-i)+i) indeed is 
empty. 

Remark: Note that variations in the energy are only possible when a particle of type 1 enters 
or leaves A or when the number of active bonds changes as a consequence of the motion of a 
particle inside A. This implies that, for all rj G Vb, ^{f],ri) for some standard configuration 
fj E 'V*j''{e*-i)+i can take only a discrete set of values. In region RB, ${f},ii) = F* can only 
happen when <?(^, -q) = H{f]) + U + 3Ai and A2 = [/ 3Ai. 

Consider the set of 2-tiled configurations consisting of a single cluster that can be reached 
with the moves considered at the various stages of the previous analysis. Clearly, this set contains 
9{v,^*,i*(i*-i)+i)- Let 77' belong to this set, and let w : — )■ 77'. Assume that max^g;^_ff(^) = 
U -|-3Ai and let supp(a;) = U5£(^supp(^). Then, from the previous analysis, it follows that there 
is a path oj' : fj ^ t]' such that max^ G uj'H{^) = H{f] + 3U and supp(cj') C supp(w). In words, if 
a configuration in T>b can be reached within energy barrier U + 3Ai, then it can also be reached 
within energy barrier 3U, irrespective of the distance from the boundary of A. ■ 

Remark 5.10 Let r] be a 2-tiled conconfiguration of minimal energy with a fixed number of 

particles of type 2, and let R{r]) be the rectangle circumscribing its dual support. Then the 2-tiled 
configuration r]' with dual support equal to R{ri) can be obtained by iteratively bringing a particle 
of type 2 to a good dual corner and saturating it with a particle of type 1 within energy barrier 
A2. Clearly, H{r}') < H{ri). We say that r( is obtained by filling the rectangle circumscribing the 
dual support of rj. 
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5.2.2 Existence of a good site 

Lemma 5.11 For all rj e Vb, there exists an x G F{'ri) such that <?((?7, x), ffl) < F*. 

Proof. The dual support of r) e Vb is either a square of side length £* or a rectangle of side 
lengths £* + — 1. Let a; be a site in a good dual corner of fj. After a particle of type 2 has 
reached site x {H{{ri, x)), T* — 3/7), it is possible to saturate this particle with an extra particle of 
type 1 within energy barrier Ai, reaching the configuration rj' with energy H{ri') = F* — A2 — e. 
Note that the dual support of 77' has the same circumscribing rectangle as ^. If the rectangle 
circumscribing the dual support of 77' is a square, then by filling this rectangle we obtain a dual 
2-tiled square of side length £*. By Remark 2.1 this is enough. If, on the other hand, the 
rectangle circumscribing the dual support of if has side length £* + I, £* — 1, then by filling this 
rectangle we obtain a dual 2-tiled rectangle of side lengths £* + 1,£* — 1. From this point on, it 
is possible to argue as in the final part of the proof of Lemma 5.7. ■ 

5.2.3 Identification of T' and Ci^d 
Lemma 5.12 V ^Vb and C^^ ^ V)^"^ . 

Proof. Same as the proof of Lemma 5.8. I 



5.3 Region RC: proof of Theorem 1.10 

Let Vc be the set of 2-tiled configurations with £*{£* — 1) + 1 particles of type 2 whose dual 
tile support is a monotone polyomino and whose circumscribing rectangle has perimeter M* (see 
Fig. 27). 



Figure 27: An example of configuration in Vc for £* = 7. 

5.3.1 Identification of g{{f^},V^,j*(^t*-i)+i) and V*,f*(^*_i)+i) 
Lemma 5.13 .g(n, K,^*(£*_i)+i) = 5(0, K,<>*(<>*_i)+i) =X>c- 

Proof. First observe that the sets V*,(?*(j?*-i)+i and "Dc coincide by Lemma 3.1. Therefore it 
remains to prove that for any configuration rj e Vc there exist a standard configuration fj and 
a path u : fj ^ rj such that max^g^j -ff (0 < or. equivalently, a path uj' : rj ^ fj such that 
max^e^' H(i) < T*. 

1. We know from Remark 4.6 that, starting from a 2-tiled configuration rj' E V^,£*(£*-i)+i 
whose dual tile support has a circumscribing rectangle of side lengths L, I with L > I and 
Lxl > — 1) + 1, it is possible to reach below energy level F* all configurations in V^,£*(£*-i)+i 



38 



whose dual support has the same circumscribing rectangle (provided that the cluster is sufficiently 
far from i9~A). We next show that, from 77', it is also possible to reach a configuration ry" £ 
V^,^£»(£._i)_i_i whose dual support has a circumscribing rectangle with side lengths L + 1,1 ~ 1 
(whenever {L + I) x {I - I) > £*{£* - 1) + 1). 

2. From 77' it is possible to reach below energy level T* a configuration 77 whose tile support, in 
dual coordinates, is a rectangle of side lengths L—1, 1 plus a bar of length k — £* {£* —1)+1—{L—1)1 
on top of the longest side of the rectangle. There are two cases. Either k > I or k < I. li k > I, 
then (L + !)(/ — 1) < £*{£* — 1) + 1. If A: < then we will show that it is possible to obtain 
within energy barrier A2 the configuration rj" whose dual tile support is obtained from the dual 
tile support of ij by moving the bar of length k from the top of the rectangle to one of its sides 
as follows. 

3. Suppose we want to move the bar onto the East side of the rectangle. Let a particle of type 
1 enter A and reach the site at dual distance 1 in the East direction from the Southern-most 
particle of type 1 on the East side of the rectangle (configuration fj') in order to create a good 
dual corner (see Fig. 28). Note that Hij)') — H{fi) + Ai. From 77', using the mechanism described 
in Section 4.1, we can iteratively move all the 2-tiles originally on the top 12-bar to the East 
side of the rectangle within energy barrier 3U. The free particle of type 1 that is left afterwards 
is removed from A. Therefore the task can be achieved within energy barrier 3U + Ai. Note that 
it is sufficient that only the North side and the East side of the dual rectangle circumscribing 
the cluster of f/ are far from A. 



'm'm'm'j ^ ^ ^ 




F — HI Si — 7* n II II ■ 1 

1- ■ 1 





Figure 28: The top bar can be moved to the East side of the rectangle within a energy barrier 
3U + Ai. 

4. It is clear that this argument is sufficient to show that from a 2-tiled configuration in 
V*^f*(f*_i)+i with a dual support consisting of a rectangle of side lengths L, (Z — 1) plus a bar 
attached to one of the shortest sides far from d^A (note that at least one of the shortest side 
is far from A) it is possible to return to a configuration of minimal energy with a dual tile 
support an L X Z rectangle and, eventually, to some standard configuration fj strictly below energy 
level r*. This implies (see Lemma 3.5) that '^(□,7;) < F* for all 77 e Dp whose dual support is 
far from A. To complete the proof we will have to consider those configurations in Vp with 
a dual tile support consisting of an L x Z rectangle that is close to d~A and see that also for 
these configurations it is possible to reach below energy level F* some standard configuration in 
Vj,_^*(^*_i)_(.i. A particle is said to be close to d^A if it is adjacent to a site in the region of A 
where interaction between particles is not possible. Dimers, 2-tiles, 12-bars and cluster are said 
to be close to d^A if they contain at least one particle that is close to d~A. 

5. We will show how it is possible to move those dimers that are close to d^A within energy 
barrier 3C/ + Ai using a modification of the argument presented in Section 4.1. For this purpose, 
we refer to Fig. 29. Let p and qi be the particle of type 2, respectively, type 1 of the dimer that is 
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(a) (b) (c) (d) 

Figure 29: Dimers that are close to A can be moved within a 3U + Ai energy barrier. A 
particle that is "beyond" the black linec an not have active bonds. 

encircled in configuration rj in Fig. 29(a). and let 52 be the particle of type 1 adjacent to p in the 
North-East direction. We will construct a path w that moves the dimer to a different 12-bar, as 
follows. Move qi one step South- East (AiJ(cj) = U), and p one step North-East {AH{uj) = 3U). 
Then move q2 one step South- East (AiJ(aj) — 3U), and let a new particle of type 1 enter A 
{AH{uj) — 3U + Ai) and reach the site originally occupied by (72 {AH{uj) = 2U + Ai; see 
Fig. 29(b)). Afterwards, move qi one step North-East {AH{uj) = 3U + Ai), p one step South- 
East (AiJ(cj) = 3C/ + Ai), and one step North-East (AiJ(w) = 2C/ + Ai; see Fig. 29(c)). The 
same procedure described so far can be repeated (now it is not necessary to let a new particle of 
type 1 enter A, since p is adjacent to two corner particles of type 1), to reach the configuration 
represented in Fig. 29(d) (as in the case of Section 4.1), which is the key configuration to see 
that the dimer can moved to a corner without further increasing the energy difference with the 
original configuration. 

6. Let rj € T>c consist of a cluster close to d~A, and let L, I be the side lengths of the rectangle 
circumscribing its dual tile support. Using the mechanism described above, we see that also from 
r] it is possible to reach a configuration 77' with a dual support consisting of a rectangle of side 
lengths £ — 1, (/) plus a 12-bar (possibly still close to d^A) attached to one of its sides. 

• li L = I = , then we can reach a standard configuration and we are done. 

• li L = + 1 and I — I* ~1 (and hence L~ I = 2), then the same procedure can be used to 
reach a configuration with a protuberance that can be easily moved below energy level F* 
to obtain a standard configuration (detach and remove the two particles of type 1, detach 
and move the particle of type 2, and saturate the particle of type 2 with two new particles 
of type 1). 

• If L — Z > 2, then it is possible to proceed as follows. Move the short external 12-bar that is 
far from 9~A onto the longest side of the dual rectangle far from d^A, as described above 
within energy barrier 3U + Ai, to obtain a configuration 77' such that H{r]') — H(ri) and 
with a dual tile support that has a circumscribing rectangle of side lengths L — 1,1 + 1. 
From rj' it is possible to iterate the above procedure until a configuration with dual support 
with a circumscribing rectangle of side lengths £* + — 1 is reached. But this case has 
already been treated. 



5.3.2 Existence of a good site 

Lemma 5.14 For all ij G Vc, there exists an x G F{fi) such that ^((t), x), ffl) < F*. 
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Proof. If 7) e Vb, then the claim follows from the proof of Lemma 5.11. We will show that, for 
all fi whose dual support has a circumscribing rectangle with side lengths L,l, there is a site x 
such that from configuration (77, x) there is a path w to a configuration rj' whose dual support 
has a circumscribing rectangle of side lengths L — 1, / + 1 such that H(fj') < H(fj) and H{^) < F* 
for all ^ G cj. The procedures that we present only require that two sides of the circumscribing 
rectangle of the dual support of t] are far from d^A, which is always the case. Without loss of 
generality we may assume that these two sides are the North side and the East side. 



Figure 30: The configuration obtained from 77 by first filling the rectangle circumscribing its dual 
support and then adding a 2-tiled on the North side; circles represent particles added to fj. 



• LI >£*{£* — 1) + 1. In this case, configuration rj contains at least one good dual corner. 
Let X be one of these dual corners. From 77, first saturate the particle of type 2 at x 
with a new particle of type 2, to obtain configuration 77'. Let 77" be the configuration 
obtained from 77' by filling i?(77) (note that 77' = 77" is possible). Clearly, H{rj") < H{t)) 
and ^(77,7;") < F*. Let 77 be the configuration obtained from rj" within energy barrier A2 
by adding a 2-tile on the North side of the 2-tiled dual rectangle with side lengths L,L 
H{fj) < F* — 4J7 + 2Ai. As in the proof of Lemma 5.13, all the tiles of the Eastern-most 
12-bar of 77 can be itereatively moved to a corner on the North side of the cluster within 
energy barrier iU (see Fig. 30), to obtain a 2-tiled configuration fj' whose dual support has 
a circumscribing rectangle with side lengths L — 1, Z + 1 and such that 77' = H{ri"). Hence 
^((7), x),f]')leT* -4U + 2Ai +3U <T* as soon as Ai < i, which is satisfied in RF. 

• LI > {£* ~ 1) + 1. In this case the dual support of 77 is already a rectangle with side 
lengths L, I. Let x be the central site of a tile adjacent to the North side of the rectangle. 
We have H{{fi,x)) = F"^ — 2U. Let 77 be the configuration obtained within energy barrier 
Ai by saturating the particle of type 2 at x. H(fj) = F* — iU + 2Ai. From 77 proceed as 
in the previous case. 



5.3.3 Identification of V and C^^ 
Lemma 5.15 V ^Vc and C^^ ^ V^'^'^. 

Proof. Same as the proof of Lemma 5.8. 
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A Computation of A^* 



Kurz [16] shows how to construct all polyominoes of minimal perimeter with fixed area, and gives 
an expression for their number in terms of generating functions. Two basic generating functions 



oo k ^ 



a{x) 



oo 



(A.l) 



are used to define two composite generating functions 

rix) = i [a{x)^ + 3a(x2)2] , q{x) = i [o(x)^ + Saix'^)'^ + 2sixfa{x^) + 2a{x^)] , (A.2) 

whose coefficients r^, qt of x'' count the polyominoes as follows. The number of polyominoes of 
minimal perimeter with area n equals 



1 



e(n) = < 



1 



Qs+l-t + Z^c=l 



[Vs+l-fJ 



if n = s^, 

if n = + i with <t < s, 
ii n = s"^ + s + t, 

i{n = s"^ + s + t with < i < s, 



(A.3) 



where s = [\/n\ . 

We need to count the number of polyominoes of minimal perimeter with area n = £*{£* — + l 
for £* > 4, i.e., we are only interested in n of the form s'^+s+t with s = £*—l and t = 1. Kurz [16] 
counts polyominoes modulo translations, rotations and reflections. We need the number modulo 
translations only. Therefore we must put in correction factors: 4 for the rotations and 2 for the 
reflections. 

In region RC we retain all c-terms. In region RB we only retain the term with c = 1. Indeed, 
qe*-i is the number of polyominoes of minimal perimeter when the circumscribing rectangle is 
a square of side length £*, and r^*_c2_i is the number of polyominoes of minimal perimeter 
when the circumscribing rectangle has side lengths £* + !,£* — 1. Thus, modulo rotations and 
reflections, we have 



8[qe*-i +n*-i-i] 



N* = 



Ec=i 



in region RB, 
in region RC. 



(A.4) 



B Clarification of some statements in [14] 
B.l Proof of Lemma 1.18 in [14] 

The following statement was used in the proof of Lemma 2.2 in [14]. 
Lemma B.l C^^ is a minimal gate. 

Proof. Let A ^ {i] e C^^: 3uj e n{r]) : uj n CJ;^ ^ {77}}, and let A = C^^\yi. In words, for each 
T] & A there is a path in (□ — ^ ffl)opt entering G{0, ffl) via 77 and reaching ffl without hitting C^^ 
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again, while if a path in (□ — )■ ffl)opt enters C^^ via a configuration in A, then it must go back 
to before reaching ffl. We will show that A is empty. The proof is by contradiction. 

Assume that r/ G A and let oj G fl{r]). Let C be the last configuration in C^^ visited by co before 
reaching ffl. Then there is an uj' E (O ^ ffl)opt entering G{0, El) via The path uj" obtained 
by joining the part of oj' from □ to ^ and the part of u) from (J to ffl belongs to (□ — )■ BEI)opt, and 
7] ^ u". Therefore r] is unessential and, by Theorem 5.1 in [18] (Lemma 1.3 in this paper), it 
does not belong to Q{0, ffl). This contradicts the assumption r] C C^^ C Q{0, ffl). ■ 

Lemma 1.18 in [14] (H3a), (H3-c) and Definition 1.5(a) imply that for every r] G C*^^ all paths 
in (ry — ^ □)opt pass through C^^. 

Proof. The proof is by contradicion. Let C*^^ 3 rj = {tItX), and assume that there is a path 
wi : □ — >■ r/ that does not visit C^^ and such that H{a) < T* for all a G uii. By the definition of 
C*t(., there exists a configuration C^^ 5 = {fj, z) such that C is obtained from r] by moving the 
particle of type 2 from x to z. Consequently, there exists a path L02 from rj to C consisting of a 
sequence of configurations of the type (r},yi) with yi G A for all i. Note that, since H{Q = T*, 
the particle of type 2 at site 2: in 1^ has no active bond (it is in d~A) and all configurations in 
UJ2 have the same number of particles of both types, we have H(a) < F* for all cr in ^2- Since C 
belongs to the minimal gate C^^, there is a path W3 : C — ^ ES such that W3 fl C^^ = {(} and such 
that H{a) < T* for all a Gus. 

Now consider the following two cases. 

(1) ri G Us. Let UJ4 be the part of W3 from r] to ffl. Then the path obtained by joining ui and 
C04 is a path in (□ ffl)opt that does not visit C^^. This contradicts the definition of C^^. 

(2) rj ^ UJ3. Let uj = uji + u>2 + wa- Then, by construction, w e (□ ^ ffl)opt- Let tt = {fj, w) 
be the configuration in oj visited just before C,. By definition, n G P. By (the first part 
of) (H3-a), TT consists of a single droplet, and hence w ^ d~A. Therefore C is obtained by 
"breaking" a droplet of a configuration in V and not by adding a particle of type 2 in d~A 
to a configuration in V. This contradicts (the second part of) (H3-a). 



B.2 Proof of Lemma 2.2 in [14] 

In the proof of Lemma 2.2 in [14], the sentence 

"Denote by 0(?7) the set of all optimal paths from □ to ffl that enter Q{D,S) via r] (note that 
this set is non-empty because C^^ is a minimal gate by Definition 1.4(a)). By Definition 1.4(b), 
cOi G ^{r]) visits fj before 77 for alH € 1, . . . , |ri(?7)|." 

must be replaced by 

"Denote by fl{r]) the set of all optimal paths from □ to ffl such that ujDC^^ = {rj} (note that this 
set is non-empty because C^^ is a minimal gate). By Definition 1.4(b) and by the more precise 
form of (H3-a), uii G ^{r]) visits f] before 77 for alH e 1, . . . , |f2(ry)|." 

B.3 Proof of Theorem 1.8 in [14] 

In Step 3, the definition of CS++(?7) should read: CS++(j?) = d+GS+ D A". 
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In Lemma 2.11, the estimate of Qi in formula (2.48) should be replaced by 

ei = [1 + 0(1)] ^ CAP^^ {d+A, CSiv)) , (B.l) 

where V is the set of configurations in P whose support has lattice distance at least two from 
d~A. This allow us to always construct a lattice path connecting every two points adjacent to 
supp(i7) and avoiding d~^A, as in the argument used to derive bounds for formula (2.56). The 
value of Oi is derived by restricting the sum in (2.50) to fj £ V . This still produces a lower 
bound for 6, because V C V. V should be changed to V accordingly in formulas (2.53) in 
Step 3 and (2.61) in Step 4. The bounds 0i and 62 will still merge asymptotically, as shown in 
Step 4, since V - iV*|A| - "P' as A ^ Z^. 

In Step 3, where Q2 is computed, the sentence 

"The only transitions in X* between C++ and U C++] are those where the free particle 

moves from distance 2 to distance 1 of the protocritical droplet." 

should be replaced by 

"The only transitions in X* between C++ and U C++] are those where the free particle 

enters CS++ (??)." 
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